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2004-2008  FINAL  PROGRESS  REPORT 

This  report  presents  the  specific  aims  and  accomplishments  of  our  breast  cancer  research  project 
during  the  years  of  funding  sponsored  by  the  US  Department  of  the  Army.  It  covers  our  activities 
from  May  1 , 2004  to  April  30,  2008. 

1.  Introduction 

The  overall  goal  of  this  research  project  is  to  apply  the  multiple  monitoring  techniques,  i.e.  Near 
infrared  spectroscopy  (NIRS),  FOXY  oxygen  sensor  and  19F  MR  EP  imaging  of  Hexafluorobenzene 
(HFB),  to  prove  the  following  hypotheses:  combination  of  hyperbaric  oxygen  (HBO)  intervention 
can  significantly  improve  breast  tumor  oxygenation,  and  that  tumor  oxygenation  remains  elevated 
for  a  substantial  period  of  time  even  after  HBO  exposure,  which  may  be  a  novel  approach  to 
enhance  radiosensitivity  or  chemotherapy.  If  our  hypotheses  are  proven  to  be  true,  this  study  can 
lead  to  an  optimal  intervention  plan  to  improve  tumor  oxygenation  and  to  determine  an  optimal  time 
interval  after  HBO  decompression  for  radiotherapy.  Such  a  novel  approach  will  largely  enhance  the 
efficiency  of  non-surgical  therapies  for  breast  tumor  treatment  and  provide  a  novel  prognostic  tool 
for  clinical  practice.  This  study  will  also  provide  a  better  understanding  of  tumor  vasculature  and 
tissue  oxygen  dynamics  and  spatial  heterogeneity  under  HBO  exposure. 

The  project  has  three  specific  aims: 

Aim  1:  to  determine  the  absolute  changes  of  oxygenated  hemoglobin  concentration,  [Hb02], 
and  hemoglobin  oxygen  saturation,  S02,  in  solid  breast  tumors  from  the  NIRS  measurements. 

Aim  2:  to  investigate  vascular  oxygenation  and  tissue  oxygen  tension  of  breast  tumors  under 
continuous  normobaric  and  hyperbaric  oxygen  exposures  with  several  gas  interventions,  using 
both  a  single-channel  NIRS  system  and  3-channel  FOXY  p02  system  simultaneously. 

Aim  3:  to  investigate  (1)  whether  HBO  could  enhance  the  therapeutic  efficiency  of  malignancy 
when  used  as  a  chemotherapeutic  adjuvant  of  an  anticancer  agent  (doxorubicin,  DOX)  in 
mammary  carcinomas  of  rat  model,  and  (2)  the  feasibility  of  NIRS  to  monitor  tumor 
hemodynamic  changes  resulting  from  the  therapeutic  effect  of  DOX  on  vasculature. 

This  report  will  summarize  the  work  we  performed  for  the  three  specific  aims  and  corresponding 
achievements,  as  given  below. 


2.  Body  of  the  Report 

The  PI  and  her  mentor,  Ms.  Mengna  Xia  and  Dr.  Liu,  made  significant  efforts  from  2004-2008  to 
accomplish  the  proposed  aims,  resulting  in  multiple  published  journal  and  conference  papers. 
This  report  is  organized  in  such  a  way  that  we  will  combine  our  published  papers  and  Dr.  Xia’s 
dissertation  to  give  detailed  descriptions  of  the  research  studies,  while  the  papers  and  chapters 
in  her  dissertation  are  organized  as  Appendixes  1  to  4  for  referral,  as  listed  below: 

Appendix  1:  Mengna  Xia.  Vikram  Kodibagkar,  Hanli  Liu,  and  Ralph  Mason,  “Tumor 
oxygen  dynamics  measured  simultaneously  by  near-infrared  spectroscopy 
and  19F  magnetic  resonance  imaging  in  rats,”  Physics  in  Medicine  and 
Biology  51 . 45-60  (2006). 

Appendix  2:  Yueqing  Gu,  Wei  R.  Chen,  Mengna  Xia.  Sang  W.  Jeong,  and  Hanli  Liu, 
“Effect  of  Photothermal  therapy  on  breast  tumor  vascular  contents:  non- 


3 


invasive  monitoring  by  near  infrared  spectroscopy,”  Photochemistry  and 
Photobiology  81  (4).  1002-1009  (2005). 


Appendix  3:  Jae  G.  Kim,  Menqna  Xia,  and  Hanli  Liu,  “Hemoglobin  extinction  coefficients: 

importance  of  correct  value  for  near-infrared  spectroscopy,”  IEEE  - 
Engineering  in  Medicine  and  Biology  magazine ,  24(2)  118-121  (2005). 

Appendix  4:  Menqna  Xia.  Hanli  Liu,  Ya  Ren,  Ralph  Mason,  Benjamin  Levine,  “Simultaneous 
monitoring  of  tumor  vascular  and  tissue  oxygen  tension  under  hyperbaric 
oxygen  exposure,”  presented  at  Biomedical  Optics  Topical  Meeting ,  Fort 
Lauderdale,  Florida,  March  19-23,  2006. 


2.1  Report  for  Aim  1 

Aim  1:  to  determine  the  absolute  changes  of  oxygenated  hemoglobin  concentration,  [Hb02], 
and  hemoglobin  oxygen  saturation,  S02,  in  solid  breast  tumors  from  the  NIRS  measurements. 

For  Aim  1,  we  have  implemented  quantification  algorithms,  examined  their  corresponding 
errors,  and  conducted  animal  measurements  to  obtain  absolute  changes  of  oxygenated 
hemoglobin  concentration,  [Hb02],  and  hemoglobin  oxygen  saturation,  S02,  in  solid  breast 
tumors  based  on  the  NIRS  measurements.  The  detailed  descriptions  can  be  found  in 

Appendixes  1  to  3. 


2.2  Report  for  Aim  2 

Aim  2:  to  investigate  vascular  oxygenation  and  tissue  oxygen  tension  of  breast  tumors  under 
continuous  normobaric  and  hyperbaric  oxygen  exposures  with  several  gas  interventions,  using 
both  a  single-channel  NIRS  system  and  3-channel  FOXY  p02  system  simultaneously. 

For  Aim  2,  we  have  reported  our  experimental  results  in  Appendixes  1  and  4  for  investigation 
of  vascular  and  tissue  oxygenation  of  rat  breast  tumors  under  continuous  normobaric  (Appendix 
1)  and  hyperbaric  (Appendix  4)  oxygen  exposures  with  oxygen  interventions,  using  single¬ 
change  NIRS  and  multi-channel  p02  sensors  simultaneously.  The  details  can  be  found  in  those 
two  appendixes. 


2.3  Report  for  Aim  3 

Aim  3:  to  investigate  (1)  whether  HBO  could  enhance  the  therapeutic  efficiency  of  malignancy 
when  used  as  a  chemotherapeutic  adjuvant  of  an  anticancer  agent  (doxorubicin,  DOX)  in 
mammary  carcinomas  of  rat  model,  and  (2)  the  feasibility  of  NIRS  to  monitor  tumor 
hemodynamic  changes  resulting  from  the  therapeutic  effect  of  DOX  on  vasculature. 

For  Aim  2,  the  PI,  Dr.  Mengna  Xia,  has  reported  the  detailed  results  in  her  dissertation.  Here  we 
report  the  major  findings  briefly,  as  follows: 

2.3.1  Introduction 


Doxorubicin  (DOX)  is  one  of  the  most  widely  used  broad-spectrum  anticancer  agents  [1], 
However,  its  clinical  utility  is  limited,  because  this  agent  produces  a  chronic  and  dose-related 
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cardiomyopathy  as  its  principal  side  effect.  Therefore,  it  is  desirable  to  achieve  better 
chemotherapeutic  effect  with  lower  dosage  of  the  agent.  It  is  well  accepted  that  hypoxic  tumor  is 
resistant  to  radiotherapy  and  some  chemotherapy  agent  [2, 3, 4, 5],  To  overcome  hypoxia,  a 
variety  of  approaches  have  focused  on  improve  oxygen  delivery  via  oxygen-enriched  gases  or 
blood  substitutes  [6,7,8,9,10].  Hyperbaric  oxygen  was  believed  to  improve  tissue  oxygenation 
greater  than  normabric  oxygen  because  it  increased  oxygen  tension  and  oxygen  delivery  to 
tissue  independent  of  hemoglobin[1 1],  HBO,  as  a  chemotherapy  adjuvant  in  tumor  treatment 
rather  than  stand  alone  treatment,  is  believed  to  increases  cellular  uptake  of  some 
chemotherapy  agents  and  the  susceptibility  of  cells  to  these  agents.  It  has  been  demonstrated 
that  HBO  can  increase  the  susceptibility  of  malignant  cells  to  destruction  with  taxol  [12], 
doxorubicin  [13,  14]  and  5-FU  [15,16], 

The  influence  of  tumor  oxygenation  on  treatment  outcome  has  stimulated  various  techniques  to 
monitor  or  estimate  tumor  oxygenation.  These  include  microelectrodes,  optical  reflectance, 
electron  paramagnetic  resonance  (EPR),  magnetic  resonance  imaging  (MRI)  and  nuclear 
medicine  approaches,  as  reviewed  previously  [17].  As  each  approach  has  their  own  strength, 
some  are  highly  invasive.  Since  its  introduction  in  1970s  [18],  Near  Infrared  Spectroscopy  has 
been  increasingly  applied  to  study  tissue  oxygenation  status  non-invasively.  Near  infrared  light 
can  easily  penetrate  biological  tissue,  and  allow  for  detection  of  specific  light-absorbing 
chromophores  in  human  in  vivo,  such  as  oxygenated  and  deoxygenated  hemoglobin,  water  and 
lipid  [19].  it  has  been  used  extensively  for  quantitative  measurements  of  cerebral  oxygenation 
[20]  and  blood  oxygenation  in  muscles  in  vivo  [21],  and  more  recently,  tumor  vascular 
oxygenation  with  respect  to  interventions  [22].  NIRS  currently  lack  of  spatial  resolution,  and 
thus,  the  utility  of  global  measurement  require  validation,  given  the  well-documented 
heterogeneity  of  tumor.  In  this  regard,  Xia  et  al  [22]  compared  the  spatially  averaged 
measurement  of  relative  tumor  oxygen  saturation  (s02)  using  NIRS  with  the  local  p02  measured 
by  MRI.  The  sensitivity  and  specificity  analysis  suggests  that  NIRS  may  identify  clinically 
relevant  hypoxia,  even  when  its  spatial  extent  is  below  the  resolution  limit  of  the  NIRS 
technique.  Kim  et  al  demonstrated  that  NIRS  may  be  used  as  an  effective  tool  to  monitor  tumor 
hemodynamic  change  induced  by  some  vascular  disrupting  agent  [23,24], 

These  studies  were  designed  to  investigate  1)  whether  HBO  could  enhance  the  therapeutic 
efficiency  of  malignancy  when  used  as  a  chemotherapeutic  adjuvant  of  doxorubicin  in  mammary 
carcinomas  of  rat  model,  and  2)  the  feasibility  of  NIRS  to  monitor  tumor  hemodynamic  changes 
resulting  from  the  therapeutic  effect  of  DOX  on  vasculature. 


2.3.2  Materials  and  Methods 


2.3.2a  Animal  and  Tumor  models 


Healthy  female  Fischer  344  rats  aged  5  to  6  months  were  obtained  from  Harlan  Sprague- 
Dawley  (Indianapolis,  IN).  They  were  housed  for  at  least  a  week  to  acclimatize  and  for 
monitoring  of  health  before  inclusion  in  the  study.  Mammary  carcinomas  13762NF  were 
implanted  in  the  dorsum  of  female  Fischer  344  rats  weighing  ~200g.  Tumor  volume  was 
estimated  by  the  formula  of  6/pi  *(LxWxH),  Tumor  diameter  was  measured  in  orthogonal  axes 
(L,W,  H).  Tumor  size  and  body  weight  were  monitored  every  other  day  after  therapeutic 
interventions.  All  animal  protocols  were  approved  by  Institutional  Animal  Care  and  Use 
Committee  at  the  University  of  Texas  Southwestern  Medical  Center  and  University  of  Texas  at 
Arlington. 
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2.3.2b  Drug  preparation  and  dose 


Doxorubicin  Hydrochloride  was  purchased  from  Sigma  Aldrich,  Inc.  It  was  made  into  a  solution 
by  dissolving  with  saline.  A  single  dose  of  DOX  (2mg/kg  body  weight)  was  administrated  by  tail 
vein.  The  dose  of  doxorubicin  given  was  the  usual  chemotherapeutic  dose,  which  would  not 
cause  cardio-toxicity  [13]. 

2.3.2c  Experimental  procedure 

Following  tumor  establishment  (~1  cm  diameter),  rats  were  randomly  assigned  to  one  of  three 
groups  according  to  different  therapeutic  strategies:  a)  DOX  (n=5),  b)  HBO  +  DOX  (n=5).  C) 
control  group  (n=2)  with  saline  injection.  Rats  were  anesthetized  with  the  mixture  of  ketamine 
hydrochloride  (0.15  ml;  lOOmg/ml;  Aveco,  Fort  Dodge,  IA)  and  xylazine  via  i.p.  After  tumor  hair 
was  shaved  to  allow  better  optical  contact  for  NIR  light  transmission,  the  rat  was  placed  on  its 
side  in  the  hyperbaric  chamber.  Probes  of  SSDRS  were  fixed  securely  on  the  tumor  of  rat,  and 
then  SSDRS  monitored  tumor  vascular  oxygenation  while  the  rat  was  subjected  to  therapeutic 
interventions. 

DOX  group  was  given  0.4  mg/ml  doxorubicin  solution  intravenously  after  the  respiratory 
intervention  of  air  -02  -  air,  and  then  were  exposed  to  air  -  02  -  air  after  DOX  injection.  HBO  + 
DOX  group  was  exposed  to  gas  intervention  in  a  sequence  of  air-02  -  HBO  (30  min)  prior  to 
intravenously  administration  of  0.4  mg/ml  doxorubicin  solution,  and  then  exposed  to  air  -  02  - 
air.  Control  group  was  subjected  to  the  respiratory  intervention  of  air  -02  -  air,  saline  solution 
intravenous  injection  and  then  were  exposed  to  air  -  02  -  air.  To  control  the  timing  of  DOX 
injection,  an  IV  butterfly  catheter  was  inserted  into  rat  tail  vein  and  fixed  securely  with  tape 
before  the  rat  was  put  into  the  chamber.  The  syringe  filled  with  DOX  solution  was  connected  to 
the  catheter  with  heparin  only  before  DOX  injection,  in  order  to  avoid  the  possible  precipitation 
caused  by  the  incompatibility  of  DOX  and  heparin  [25]. 

2. 3. 2d  Spectrometer  for  monitoring  the  disturbance  of  DOX  on  the  absorption  of  tissue 
phantom 

Because  the  reddish  color  of  doxorubicin,  it  is  reasonable  to  consider  the  bolus  injection  of 
doxorubicin  may  change  tissue  absorption.  A  UV/VIS  spectrometer  Lambda  20  (PerkinElmer 
Inc.,  Waltham,  MA)  was  used  to  detect  the  absorption  change  resulting  from  the  addition  of 
DOX  into  tissue  phantom,  in  order  to  verify  the  acute  effect  of  DOX  bolus  injection  on  tissue 
absorption.  A  tissue  phantom  is  composed  of  100  I  sheep  blood  mixed  with  phosphate  buffer 
solution  with  total  volume  of  3.5  ml  and  total  hemoglobin  concentration  of  7.1  g/L,  which  is  in  the 
range  of  total  hemoglobin  concentration  in  tissue  [26].  The  Hb  concentration  of  tissue  phantom 
was  measured  by  Co-oximeter  (Instrumentation  Lab,  Ramsey,  MN).  Since  animal’s  total  blood 
volume  is  10%  of  its  body  weight,  total  blood  volume  =  0.2  Kg  x  10%  =  20  ml,  with  the  known 

body  weight  of  ~  0.2  Kg  of  rats.  0.17  ml  DOX  solution  with  the  concentration  of  0.4  mg/ml  was 

added  into  the  tissue  phantom.  Accordingly,  the  volume  ratio  of  DOX  to  tissue  phantom  is 
proportional  to  the  ratio  of  1  ml  0.4  mg/ml  DOX  to  20  ml  total  blood  volume  in  rats.  According  to 
the  Lambert  Beer  Law,  the  absorption  in  750nm  and  830nm  are 

A 750  (0)  =  sHb 750  [Hb]  +  sHb0i  750  [Hb02  ]  (2.1) 

A750  (1)  =  sHb750  [Hb]  +  sHWi  750  [Hb02  ]  +  sDOX150  [DOX]  (2.2) 

^83°(0)  =  sH™[Hb]  +  sHh0™[Hb02]  (2.3) 


6 


(2.4) 


a  830  /i  \  830  r-  yti  830  r  yti  ~i  830  r  t/--i 

A  (1)  =  s  Hb  \Hb ]  +  £  HbQ2  \Hb02  ]  +  £  DOx  \DOX\ 

where  A750  (0)  and  A750  (1)  represent  the  absorption  at  750  nm  in  tissue  phantom  without  DOX 
and  with  DOX,  respectively.  A830  (0)  and  A830  (1  )  represent  the  absorption  at  830  nm  in  tissue 
phantom  without  DOX  and  with  DOX,  respectively. 

2.3.2e  Steady-state  diffuse  reflectance  spectroscopy  (SSDRS)  for  measuring  changes  in  tumor 
vascular  oxygenation  ( AHbO ?) 

A  broadband  diffuse  reflectance  spectrometer  was  used  to  acquire  reflectance  spectra  from 
tumor  tissue.  Briefly,  continuous  wave  (CW)  light  from  a  20  W  tungsten-halogen  light  source 
(HL-2000HP,  ocean  optics,  FL)  is  coupled  into  a  2.6-mm  core  diameter  fiber  optic  bundle,  the 
distal  end  of  which  is  placed  in  physical  contact  with  the  surface  of  the  tumor.  After  being 
scattered  in  the  tumor  tissue,  the  transmitted  light  is  collected  by  a  1-mm  core  diameter 
detection  fiber,  the  end  of  which  is  coupled  to  a  hand-held  spectrometer  (USB2000,  Ocean 
optics,  FL).  The  broadband  light  diffuse  spectrometer  provides  reflectance  spectra  from  400  to 
900  nm. 


According  to  the  modified  Beer-Lambert  law,  changes  of  oxy-  and  deoxy-hemoglobin 
concentration,  A[Hb02]  and  A[Hb],  can  be  derived  from  the  measured  amplitudes  at  two 
wavelengths  (750nm  and  830nm),  by  using  extinction  coefficients  of  oxy-  and  deoxy-Hb 
published  by  Cope  [27],  as  given  in  Equations.  (2.5)  and  (2.6): 


A  [HbQ2\ 


-1.532  log(^)750  + 1.753  log(^) 830 

_ A _ A 

DPF  -d 


A  [Hb]  = 


1.758  log( ^-) 750  -  0.92  log(— )’ 330 

_ A _ A 

DPF  -d 


(2.5) 


(2.6) 


where  Ab  is  the  baseline  amplitude,  At  is  the  transient  amplitude  during  the  intervention,  and  d  is 
the  direct  source-detector  separation.  DPF  (differential  path-length  factor)  is  a  tissue-dependent 
parameter  and  defined  as  the  ratio  between  the  optical  path  length  and  the  physical  separation 
between  the  source  and  detector. 

2.3.3  Results 


2.3.3a  Disturbance  of  DOX  on  the  absorption  of  tissue  phantom 

There  is  absorption  difference  between  tissue  phantom  with  and  without  DOX,  as  shown  in 
Figure  2.1.  However,  the  differences  are  relatively  small.  Indeed,  the  relative  changes  of 

A 730  (1) _ A 730  (0) 

absorption  in  wavelength  of  750nm  and  830nm  are  - — - =4.85%  and 

A750  (  0) 
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^4830m  _  ^4830  (  0) 

- — — - =4.80%,  respectively.  In  order  to  investigate  the  disturbance  of  DOX  on 

^830  (  0) 

calculating  oxygen  saturation,  absorption  spectrum  profiles  between  700  and  900  nm  were 
normalized  by  absorption  values  at  the  wavelength  of  700  nm.  The  normalized  spectra  cover 
the  wavelengths  we  utilized  for  calculating  hemoglobin  concentration.  As  shown  in  figure  2.2, 
profiles  of  normalized  absorption  spectrum  in  NIR  range  appeared  to  be  overlaid. 


Wavelength  (nm) 


Figure  2.1  Absorption  spectra  of  tissue  phantoms  with  (yellow  curve)  or  without  1 70  I  DOX 
(blue  curve).  The  unit  for  absorption  is  arbitrary  unit. 


Wavelength  (nm) 

Figure  2.2  Normalized  absorption  spectra  of  tissue  phantoms  with  (yellow  curve)  or  without 
170  I  DOX  (blue  curve)  in  the  wavelength  range  of  700  ~  900  nm. 


2.3.3b  Changes  in  Tumor  volume  and  body  weight  during  chemotherapy 

Tumor  volume  and  body  weight  were  monitored  before  and  after  DOX  treatment  to  examine  the 
tumor  response.  Changes  in  tumor  volume  and  body  weight  were  normalized  by  the  values  at 
day  0  (before  DOX  administration). 


Day 

Figure  2.3  Normalized  (a)  tumor  volume  and  b)  body  weight  in  rats  with  saline  injection  (•), 
DOX  treatment  (  □  )  and  HBO  +  DOX  treatment  (  A  ). 


Tumors  in  the  control  group  grow  significantly  faster  than  tumors  with  treatment  in  the  other 
groups  (p<0.05),  as  shown  in  Figure  2.3a.  Tumors  with  combined  treatment  of  HBO  and  DOX 
grow  significantly  slower  than  those  with  DOX  treatment  except  for  the  first  two  days  after 
treatment  (p<0.05).  Basically,  there  is  no  significant  difference  for  tumor  volume  in  HBO  +  DOX 
and  DOX  groups  on  the  2nd  day  after  treatment.  However,  significant  differences  are  observed 
starting  from  5th  day  after  treatment  between  the  DOX  treated  and  control  group  (p<  0.05). 
Regarding  body  weight  loss,  figure  2.3b  indicated  that  rats  in  DOX  group  and  DOX  +  HBO 
group  had  significant  and  continuous  body  weight  loss  after  day  0.  Rats  in  control  group  gained 
weight  at  day  2,  and  started  to  lose  weight  after  day  4,  and  then  kept  constant  body  weight 
thereafter. 

2.3.3c  Vascular  hemodynamic  changes  of  rats  in  DOX  group 
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Values  of  A[Hb02]  of  rats  in  DOX  group  showed  increases  with  values  in  a  range  of  0.07  -0.12 
(mM/DPF),  when  the  gas  switched  from  air  to  oxygen  before  DOX  injection.  When  the  gas  was 
switched  from  air  to  oxygen  after  DOX  administration,  the  maximal  amplitude  of  A[Hb02] 
increased  with  values  in  a  range  of  0.02  -  0.04  (mM/DPF),  which  is  significantly  less  than  the 
increase  achieved  with  02  inhalation  before  DOX  administration,  as  shown  in  figure  2.4.  We 
also  noticed  the  signal  changes  during  DOX  i.v.  injection. 


Figure  2.4  Time  profile  of  A[Hb02]  of  tumors  (1  -4)  in  DOX  group  when  the  rat  was  under 
gas  intervention.  The  unit  for  A[Hb02]  is  mM/DPF.  A[Hb02]  of  the  5th  rat  was  discarded 
because  of  misplaced  gas  mask  during  the  experiment  due  to  the  movement  of  the  rat. 


to 


2.3.3d  Vascular  hemodynamic  changes  of  rats  in  POX  +HBO  group 


Time»  /mini 


Time  (min) 


Time  (min) 


Time  (min) 


Time  (min) 


Figure  2.5  Typical  time  profile  of  A[Hb02]  of  tumors  in  DOX  +  HBO  group  (n=5) 
when  the  rat  was  under  gas  intervention. 
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It  shows  a  stable  baseline  in  A[Hb02]  when  rats  were  inhaling  air.  A[Hb02]  increase  immediately 
in  the  first  few  minutes  and  more  gradually  afterwards,  as  shown  in  figure  2.5.  When  the  rats 
were  exposed  to  hyperbaric  oxygen,  A[Hb02]  has  a  further  increase  until  reaching  a  stabilized 
value.  DOX  i.v.  injection  immediately  after  HBO  caused  some  fluctuation  of  signal,  but  stabilized 
when  the  injection  is  finished.  A[Hb02]  has  a  stable  baseline  but  with  values  greater  than  air 
inhalation  before  HBO  when  the  rats  were  breathing  air  after  HBO.  Similar  to  DOX  group,  the 
maximal  amplitude  of  A[Hb02]  achieved  with  02  inhalation  prior  to  the  DOX  administration  is 
significantly  greater  than  that  achieved  with  02  inhalation  after  DOX  administration. 


2.3.3e  Vascular  hemodynamic  changes  of  rats  in  control  group 

A[Hb02]  have  stable  baseline  values  when  rats  were  breathing  air,  as  shown  in  figure  2.6. 
A[Hb02]  increased  when  the  gas  was  switching  to  oxygen,  and  decreased  when  the  gas  was 
switched  back  to  air,  as  expected.  Different  from  A[Hb02]  in  DOX  treated  group,  the  change  of 
amplitude  in  A[Hb02]  when  the  gas  switched  from  air  to  oxygen  after  saline  injection  is  similar  to 
the  change  due  to  gas  intervention  before  injection  in  both  rats. 
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Figure  2.6  Time  profile  of  A[Hb02]  of  tumors  in  control  group  (n=2)  when  rats  were 
subjected  to  gas  intervention. 


2.3.4  Discussion  and  Conclusion 


In  our  study,  we  randomized  breast  tumors  into  three  groups:  group  with  DOX  injection,  group 
with  DOX  injection  following  HBO  exposure  at  2  atm,  control  group  with  saline  injection.  The 
result  showed  that  tumors  in  HBO  +  DOX  group  grow  significantly  slower  than  those  with  DOX 
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alone.  HBO  enhanced  chemotherapeutic  response  of  mammary  carcinoma  NF13762  to  DOX  in 
vivo  reflected  by  slowing  down  the  tumor  growth  after  treatment. 

Resistance  to  chemotherapy  is  common  in  hypoxic  tumors.  HBO  may  help  overcome 
chemotherapy  resistance  by  improving  both  tumor  perfusion  and  cellular  sensitivity.  Improving 
tumor  oxygenation  and  vascularization  may  increase  drug  delivery.  This  has  been  shown 
experimentally  and  in  nude  mice  with  human  epithelial  ovarian  cancer  treated  with  cisplatin  [28], 
Reactive  oxygen  species  (ROS),  or  free  radicals  are  by-product  of  aerobic  respiration  and 
cellular  metabolism  and  induced  by  oxidative  stress  during  hypoxia  (oxygen  deficiency), 
reperfusion  or  hyperoxia  (excess  oxygen).  ROS,  at  low  levels,  assist  tumor  growth  but  become 
toxic  at  high  levels.  This  can  be  explained  by  the  “threshold  effect”  whereby  ROS  reach  a  level 
beyond  which  the  antioxidant  capacity  is  inundated,  resulting  in  irreversible  damage  and 
apoptosis  [29,30].  One  of  the  mechanisms  of  action  of  doxorubicin  is  production  of  ROS.  By 
increasing  ROS  level,  HBO  push  ROS  levels  past  the  threshold  level,  and  thus  enhanced  the 
ROS-localized  effects  of  doxorubicin  [31 ,32],  Another  mechanism  of  HBO  is  to  push  cell  to  enter 
a  proliferate  stage,  thus  sensitizing  them  to  radiotherapy  and  some  chemotherapy  by  improving 
oxygenation.  It  has  been  showed  that  HBO  enhanced  the  chemotherapeutic  effects  of 
doxorubicin  in  an  experimental  model  of  pulmonary  sarcoma  [13].  HBO  stimulated  proliferation 
of  an  MCA-2  metastatic  lung  tumor  cell  line  and  induced  cells  to  enter  the  replicating  cycle 
compared  to  cells  left  at  ambient  pressure  [13].  Other  studies  found  that  HBO  increased  the 
ratio  of  prostate  cancer  cells  in  vitro  accumulating  in  G2/M  phases  from  the  G0  arrest  phase  [33]. 
Generally  speaking,  HBO  therapy  in  combination  with  chemotherapy  may  be  justified  by  the 
following:  1)  improved  oxygenation  improves  drug  delivery  to  hypoxic  regions  in  the  tumor;  2) 
increasing  intratumoral  ROS  levels  beyond  the  threshold  may  induce  tumor  destruction;  3) 
improved  oxygenation  may  also  cause  cell  to  enter  a  proliferate  stage,  thus  sensitizing  them  to 
radiotherapy  and  some  chemotherapy;  4)  HBO  may  remove  hypoxia  stimulus  that  drives 
angiogenesis. 

Even  though  the  DOX  dosage  is  reported  to  cause  minimal  cardiotoxicity,  the  result  showed  that 
changes  in  amplitudes  of  A[Hb02]  before  DOX  administration  were  much  greater  than  the 
change  after  treatment  in  rats  of  DOX  group  and  DOX  +  HBO  group.  It  is  likely  that  the 
amplitude  difference  results  from  the  known  cardiotoxicity  reaction,  the  major  side  effect  of 
DOX.  It  has  been  suggested  that  cardiac  dysfunction  induced  by  DOX  resulted  from  the 
imbalances  of  the  circulatory  system  such  as  decreases  in  blood  pressure  or  the  direct  effects 
on  vascular  wall  [34].  During  the  course  of  i.v.  injection  of  DOX,  the  vasculature  was  exposed  to 
high  levels  of  DOX,  and  in  vitro  studies  have  suggested  that  DOX  acutely  induces  vascular 
smooth  muscle  to  release  Ca2+  from  its  intracellular  storage  site  and  causes  direct 
vasoconstrictor  [35]  and  vasodilator  effects  [36],  Furthermore,  the  combination  of  doxorubicin 
and  HBO  would  also  be  expected  to  enhance  the  agent’s  cardiotoxicity  because  of  the  toxicity 
to  cardiomyocytes  of  HBO.  Therefore,  the  clinical  addition  of  HBO  to  doxorubicin  may  not 
change  the  risk-benefit  ratio  of  the  agent.  NIRS,  in  turn,  may  provide  a  novel  approach  to 
monitor  the  cardiotoxicity  of  treatment,  which  may  leads  to  an  optimized  therapeutic  plan  to 
minimize  the  side  effect  of  treatment.  We  also  noticed  the  signal  fluctuation  in  A[Hb02]  during 
DOX  injection  in  DOX  group.  DOX  solution  is  orange-red,  so  it  is  likely  that  DOX  bolus  injection 
would  cause  the  absorption  change  of  tumor  tissue  in  NIR  range.  We  measured  and  compared 
the  spectra  of  tissue  phantoms  before  and  after  adding  DOX,  to  examine  the  disturbance  of 
DOX  injection  on  the  signal  of  A[Hb02].  There  were  absorption  differences  when  comparing 
both  absorption  spectra  (Figure  2.1).  However,  the  normalized  spectra  appeared  to  be  overlaid 
(Figure  2.2).  Therefore,  it  implied  that  DOX  would  affect  total  hemoglobin  concentration,  rather 
than  oxygen  saturation. 
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3.  Key  Research  Accomplishments  and  Reportable  Outcomes 

3.1  Ph.D.  Dissertation: 

Title:  TISSUE  AND  VASCULAR  OXYGENATION  DYNAMICS  DETERMINED  BY 
OPTICAL  APPROACHES  AND  MRI 

3.2  Publications  in  peer-reviewed  journals: 

(1)  Mengna  Xia,  Vikram  Kodibagkar,  Hanli  Liu  and  Ralph  Mason,  “Tumor  oxygen  dynamic 
measured  simultaneously  by  near  infrared  spectroscopy  and  19F  magnetic  resonance 
imaging  in  rats”,  Physics  in  Medicine  and  Biology, 51:  45-60  (2006) . 

(2)  Jae  Kim,  Mengna  Xia,  Hanli  Liu,  “Extinction  coefficients  of  hemoglobin  for  near-infrared 
spectroscopy  of  tissue”,  IEEE  Engineering  in  Medicine  and  Biology  magazine,  24:  118-121 
(2005). 

(3)  Yueqing  Gu,  Wei  R.  Chen,  Mengna  Xia,  Sang  W.  Jeong,  and  Hanli  Liu,  “Effect  Of 
photothermal  therapy  On  breast  tumor  vascular  contents:  non-invasive  monitoring  by  near 
infrared  spectroscopy”,  Photochemistry  and  Photobiology,  81:  1002-1009  (2005). 

3.3  Proceeding  papers  and  presentations: 

(1)  Mengna  Xia,  Benjamin  Levine,  Ralph  Mason,  Hanli  Liu,  “Simultaneous  monitoring  of  tumor 
vascular  oxygenation  and  tissue  oxygen  tension  under  hyperbaric  oxygen  exposure”,  in 
Biomedical  Topical  Meetings  on  CD-ROM  (The  Optical  Society  of  America,  Washington, 

DC,  2006). 

(2)  Mengna  Xia,  Vikram  Kodibagkar,  Ralph  Mason,  Benjamin  Levine,  Hanli  Liu,  “Tumor 
vascular  and  tissue  oxygenation  dynamics  under  normobaric  and  hyperbaric  oxygen 
interventions”,  presented  at  the  fourth  Era  of  Hope  meeting  for  the  Department  of  Defense 
(DOD)  Breast  Cancer  Research  Program  (BCRP)  held  on  June  8-11,  2005  in  Philadelphia, 
Pennsylvania. 

(3)  Mengna  Xia,  Ralph  Mason,  Hanli  Liu,  “A  model  of  hemodynamic  responses  of  rat  tumors  to 
hyperoxic  gas  challenge”,  Proc.  SPIE-  Optical  Tomography  and  Spectroscopy  of  Tissue  VI, 
5693:  301-307  (2005). 

4.  Conclusions 

In  this  4-year  study,  we  have  learned  that  HBO  enhances  the  therapeutic  action  of  doxorubicin 
in  our  tumor  model,  probably  by  multiple  physiological  mechanisms.  The  present  study  reveals 
that  DOX  may  be  used  in  conjunction  with  HBO  to  obtain  the  same  effect  as  higher  doxorubicin 
doses.  Meanwhile,  NIRS  may  work  as  an  attractive  approach  to  monitor  the  treatment. 

Overall,  our  study  over  the  last  few  years  has  demonstrated  that  1)  NIRS  and  MRI  and  FOXY™ 
oxygen  sensor  are  complimentary  approaches  to  monitor  tumor  oxygenation.  2)  Tumor  tissue 
oxygenation  achieved  by  hyperbaric  oxygenation  persists  over  10-20  minutes  even  after 
terminating  hyperbaric  oxygenation  intervention.  3)  Several  correlations  existed  for  both 
modalities  (Hb02  and  p02)  under  sequences  of  hyperoxic  gas  intervention  with  hyperbaric 
oxygen  exposure.  Correlation  of  tumor  vascular  oxygenation  and  tumor  tissue  p02  determined 
simultaneously  by  those  techniques  could  give  us  a  better  understanding  on  the  patho¬ 
physiology  of  tumor  and  response  to  therapeutic  interventions.  4)  HBO  enhanced  the 
therapeutic  action  of  doxorubicin  in  our  tumor  model,  probably  by  multiple  physiological 
mechanisms.  DOX  could  be  used  in  conjunction  with  HBO  to  achieve  the  same  effect  as  higher 
doxorubicin  doses.  NIRS  may  work  as  an  attractive  approach  to  monitor  radiation-toxicity  of  the 
treatment. 


14 


5.  References 


1  .  Weiss,  R.,  G.  Sarosy,  K.  Clagett-carr,  M.  Russo,  and  B.  Leyland-jones,  Anthracyline 
analogs:  the  past,  present,  and  future.  Cancer  Chemother  Pharmacol,  1986.  18:  p.  185-97. 

2.  Brown,  J.M.,  The  hypoxic  cell:  a  target  for  selective  cancer  therapy-eighteenth  Bruce  F. 

Cain  Memorial  Award  lecture.  Cancer  Res.,  1999.  59(23):  p.  5863-5870. 

3.  Hall,  E.,  Radiobiology  for  the  Radiologist.  4th  ed.  1994,  Philadelphia,  PA:  Lippincott. 

4.  Thews,  O.,  D.  Kelleher,  and  P.  Vaupel,  Erythropoietin  restores  the  anemia-induced  reduction 
in  cyclophosphamide  cytotoxicity  in  rat  tumors.  Cancer  Res,  2001.  61:  p.  1358-61. 

5.  Bush,  R.,  R.  Jenkin,  W.  Allt,  F.  Beale,  A.  Dembo,  and  J.  Pringle,  Definitive  evidence  for 
hypoxic  cells  influencing  cure  in  cancer  therapy.  Br  J  Cancer,  1978.  37:  p.  302-06. 

6.  Hall,  E.J.,  The  oxygen  effect  and  reoxygenation,  in  Radiobiology  for  the  Radiologist,  E.J. 

Hall,  Editor.  1994,  J.  B.  Lippincott:  Philadelphia,  p.  133-152. 

7.  Brady,  L.,  H.  Plenk,  J.  Hanley,  J.  Glassburn,  S.  Kramer,  and  R.  Parker,  Hyperbaric  oxygen 
therapy  for  carcinoma  of  the  cervix  stages  I  IB,  MB,  and  IVA:  results  of  a  randomized  study 
by  the  radiation  therapy  oncology  group.  Int  J  Radiat  Onclo  Biol  Phys,  1981 . 7:  p.  991-8. 

8.  Brizel,  D.,  W.  Hage,  R.  Dodge,  M.  Munley,  C.  Piantadosi,  and  M.  Dewhirst,  Hyperbaric 
Oxygen  Improves  Tumor  Radiation  Response  Significantly  More  Than 
carbogen/Nicotinarnid.  Radiation  Res.,  1997.  147:  p.  715-20. 

9.  Rockwell,  S.,  M.  Kelley,  C.  Irvin,  C.  Hughes,  E.  Porter,  H.  Yabuki,  and  J.  Fischer,  Modulation 
of  tunior  oxygenation  and  radiosensitivity  by  a  perfluorooctylbromide  emulsion.  Radiother 
Oncol,  1991. 22:  p.  92-8. 

10.  Lustig,  R.,  N.  Lowe,  C.  Rose,  J.  Haas,  S.  Krasnow,  M.  Spaulding,  and  I.  Prosnitz,  Phase 
1/11  study  kluosol  and  100%  oxygen  as  an  adjuvant  to  radiation  ill  the  treatment  of  locally 
advanced  non-small  cell  carcinoma  of  the  head  and  neck.  Int.  .1.  Radiat.  Oncal.  Biol.  Phys., 
1989.  16:  p.  1587-93. 

11.  Xia,  M.,  V.  Kodibagkar,  R.  Mason,  B.  Levine,  and  H.  Liu.  Tumor  vascular  and  tissue 
oxygenation  dynamics  under  normobaric  and  hyperbaric  oxygen  interventions,  in  fourth  Era 
of  Hope  meeting  for  the  Department  of  Defense  (DOD)  Breast  Cancer  Research  Program 
(BCRP).  2005.  Philadelphia,  PA. 

12.  Kalns,  J.,  L.  Krock,  and  E.  Piepmeier,  The  effect  of  hyperbaric  oxygen  on  growth  and 
chemosensitivity  of  metastatic  prostate  cancer.  Anticancer  Res,  1998.  18. 

13.  Petre,  P.,  F.  Baciewicz,  S.  Figan,  and  J.  Spear,  Hyperbaric  oxygen  as  a  chemotherapy 
adjuvant  in  the  treatment  of  metastatic  lung  tumors  in  a  rat  model.  J  Thorac  Cardiovasc 
Surg,  2003.  125:  p.  85-95. 

14.  Kalns,  J.,  L.  Krock,  and  E.  Piepmeier,  The  effect  of  hyperbaric  oxygen  on  growth  and 
chemosensitivity  of  metastatic  prostate  cancer.  Anticancer  Res,  1998.  18. 


15 


15.  Takiguchi,  N.,  N.  Saito,  M.  Nunomura,  K.  Kouda,  K.  Oda,  N.  Furuyama,  and  N.  Nakajima, 
Use  of  5-FU  plus  hyperbaric  oxygen  for  treating  malignant  tumorsievaluation  of  antitumor 
effect  and  measurement  of  5-FU  in  individual  organs.  Cancer  Chemother  Pharmacol,  2001 . 
47:  p.  11-14. 

16.  Stuhr,  L.,  Hyperbaric  oxygen  alone  or  combined  with  5-FU  attenuates  growth  of  DMBA- 
induced  rat  mammary  tumors.  Cancer  Lett,  2004.  210:  p.  35-40. 

17.  Zhao,  D.,  L.  Jiang,  and  R.  Mason,  Measuring  changes  in  tumor  oxygenation.  Methods 
Enzymol,  2004.  386:  p.  378-418. 

18.  Jobsis,  F.,  Noninvasive,  infrared  monitoring  of  cerebral  and  myocardial  oxygen  sufficiency 
and  circulatory  parameters.  Science,  1977.  198:  p.  1264-7. 

19.  Liu,  H.,  Y.  Song,  K.L.  Worden,  X.  Jiang,  A.  Constantinescu,  and  R.P.  Mason,  Noninvasive 
Investigation  of  Blood  Oxygenation  Dynamics  of  Tumors  by  Near-Infrared  Spectroscopy. 
Appl.  Opt.,  2000.  39(28):  p.  5231-5243. 

20.  Delpy,  D.T.  and  M.  Cope,  Quantification  in  tissue  near-infrared  spectroscopy.  Phil.  Trans. 
R.  Soc.  Lond.  B„  1997.  352:  p.  649-659. 

21.  Homma,  S.,  T.  Fukunaga,  and  A.  Kagaya,  Influence  of  adipose  tissue  thickness  on  near- 
infrared  spectroscopic  signals  in  the  measurement  of  human  muscles.  J.  Biomed.  Opt., 
1996.  1:  p.  418-424. 

22.  Xia,  M.,  V.  Kodibagkar,  H.  Liu,  and  R.P.  Mason,  Tumour  oxygen  dynamics  measured 
simultaneously  by  near-infrared  spectroscopy  and  19F  magnetic  resonance  imaging  in  rats. 
Physics  in  medicine  and  biology,  2006.  51:  p.  45-60. 

23.  Kim,  J.,  D.  Zhao,  R.  Mason,  and  H.  Liu.  Acute  effects  of  combreatastatin  A4  phosphate  on 
breast  tumor  hemodynamics  monitored  by  Near  Infrared  Spectroscopy,  in  in  Biomedical 
Optics  2006  Technical  Digest  (Optical  Society  of  America,  Washington,  DC,  2006).  2006. 

24.  Kim,  J.,  D.J.  Cuccia,  J.  Lee,  A.E.  Cerussi,  A.J.  Durkin,  and  B.J.  Tromberg.  Monitoring  of 
chemotherapy  effects  on  rat  breast  tumors:  a  correlation  between  optical  signals  and 
histology,  in  SPIE.  2007.  San  Jose,  CA. 

25.  Cohen,  M.,  A.  Johnston-Early,  M.  Hood,  M.  McKenzie,  M.  Citron,  N.  Jaffe,  and  S. 

Krasnow,  Drug  precipitation  within  i.v.  tubing:  a  potential  hazard  of  chemotherapy 
administration.  Cancer  Treat  Rep.  1985  Nov;69(1 1):  1325-6,  1985.  69:  p.  1325-6. 

26.  Esenaliev,  R.,  Y.  Petrov,  O.  Hartrumpf,  D.  Deyo,  and  D.  Prough,  Continuous,  noninvasive 
monitoring  of  total  hemoglobin  concentration  using  an  optoacoustic  technique.  Appl  Opt, 
2004.  43:  p.  3401-07. 

27.  Cope,  M.,  the  application  of  near  infrared  spectroscopy  to  non  invasive  monitoring  of 
cerebral  oxygenation  in  the  newborn  infant.  1991,  Ph.D  dissertation  in  University  of  London. 


16 


28.  Alagoz,  T.,  R.  Buller,  B.  Anderson,  K.  Terrell,  R.  Squatrito,  T.  Niemann,  D.  Tatman,  and  P. 
Jebson,  Evaluation  of  hyperbaric  oxygen  as  a  chemosensitizer  in  the  treatment  of  epithelial 
ovarian  cancer  in  xenografts  in  mice.  Cancer,  1995.  75:  p.  2313-22. 

29.  Kong,  Q.,  J.  Beel,  and  K.  Lillehei,  A  threshold  concept  for  cancer  therapy.  Med  Hpyotheses, 
2000.  55:  p.  29-35. 

30.  Harrisona,  L.,  M.  Chadhaa,  R.  Hillb,  K.  Hua,  and  D.  Shashaa,  Impact  of  tumor  hypoxia  and 
anemia  on  radiation  therapy  outcomes.  Oncologist,  2002.  7:  p.  492-508. 

31.  Kizaka-Kondoh,  S.,  M.  Inoue,  H.  Harada,  and  M.  Hiraoka,  Tumor  hypoxia:  a  target  for 
selective  cancer  therapy.  Cancer  Sci,  2003.  94:  p.  1021-28. 

32.  McMillan,  T.,  K.  Calhoun,  J.  Mader,  C.  Stiernberg,  and  S.  Rajaraman,  The  effect  of 
hyperbaric  oxygen  on  oral  mucosal  carcinoma.  Laryngoscope,  1989.  99:  p.  241-44. 

33.  Kalns,  J.  and  E.  Piepmeier,  Exposure  to  hyperbaric  oxygen  induces  cell  cycle  perturbation 
in  prostate  cancer  cells.  In  Vitro  Cell  Dev  Biol  Anim,  1999.  35:  p.  98-101. 

34.  Murata,  T.,  H.  Yamawaki,  M.  Hori,  K.  Sato,  H.  Ozaki,  and  H.  Karaki,  Chronic  vascular 
toxicity  of  doxorubicin  in  an  organ-cultured  artery.  British  J  of  Pharm,  2001. 132:  p.  1365-73. 

35.  Kanmura,  Y.,  L.  Raeymaekers,  and  R.  Casteels,  Effects  of  doxorubicin  and  ruthenium  red 
on  intracellular  Ca2+  stores  in  skiined  rabbit  mesenteric  smooth-muscle  fibres.  Cell  Calcium, 
1989.  10:  p.  433-39. 

36.  Ferrans,  V.,  J.  Clark,  J.  Zhang,  Z.  Yu,  and  E.  Herman,  Pathogenesis  and  prevention  of 
doxorubicin  cardiomyopathy.  Tsitologiia,  1997.  39:  p.  928-37. 


17 


Institute  of  Physics  Publishing 


Physics  in  Medicine  and  Biology 


Phys.  Med.  Biol.  51  (2006)  45-60 


doi:  10. 1088/003 1  -9 155/5 1/1/004 


Tumour  oxygen  dynamics  measured  simultaneously 
by  near-infrared  spectroscopy  and  19F  magnetic 
resonance  imaging  in  rats 


Mengna  Xia1,  Vikram  Kodibagkar2,  Hanli  Liu1  and  Ralph  P  Mason1,2 

1  Joint  Biomedical  Engineering  Graduate  Program,  University  of  Texas  at  Arlington/University 
of  Texas  Southwestern  Medical  Center  at  Dallas,  TX  76019,  USA 

2  Cancer  Imaging  Program,  Department  of  Radiology,  University  of  Texas  Southwestern  Medical 
Center  at  Dallas,  TX  75390,  USA 

E-mail:  Ralph.Mason@UTSouthwestern.edu 

Received  1  April  2005,  in  final  form  17  October  2005 
Published  14  December  2005 
Online  at  stacks.iop.org/PMB/51/45 

Abstract 

Simultaneous  near-infrared  spectroscopy  (NIRS)  and  magnetic  resonance 
imaging  (MRI)  were  used  to  investigate  the  correlation  between  tumour 
vascular  oxygenation  and  tissue  oxygen  tension  dynamics  in  rat  breast  1 3762NF 
tumours  with  respect  to  hyperoxic  gas  breathing.  NIRS  directly  detected 
global  variations  in  the  oxygenated  haemoglobin  concentration  (AfHbOi]) 
within  tumours  and  oxygen  tension  (pOj)  maps  were  achieved  using  19F 
MRI  of  the  reporter  molecule  hexafluorobenzene.  Multiple  correlations  were 
examined  between  rates  and  magnitudes  of  vascular  (A[HbOi])  and  tissue 
(pCF)  responses.  Significant  correlations  were  found  between  response  to 
oxygen  and  carbogen  breathing  using  either  modality.  Comparison  of  results 
for  the  two  methods  showed  a  correlation  between  the  vascular  perfusion 
rate  ratio  and  the  mean  pCF  values  (R2  >  0.7).  The  initial  rates  of  increase 
of  A[  HbCU  |  and  the  slope  of  dynamic  pCU  response,  d(p()ij/d/.  of  well- 
oxygenated  voxels  in  response  to  hyperoxic  challenge  were  also  correlated. 
These  results  demonstrate  the  feasibility  of  simultaneous  measurements  using 
NIRS  and  MRI.  As  expected,  the  rate  of  pCF  response  to  oxygen  is  primarily 
dependent  upon  the  well  perfused  rather  than  poorly  perfused  vasculature. 


1.  Introduction 

Tumour  oxygenation  has  been  widely  recognized  as  a  pivotal  factor  in  the  efficacy 
of  radiotherapy  (Hall  1994),  photodynamic  therapy  (Chapman  et  al  1991)  and  some 

*  Presented  in  part  at  the  12th  annual  meeting  of  the  International  Society  of  Magnetic  Resonance  in  Medicine, 
Kyoto,  2004. 
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chemotherapies  (Brown  1999),  and  patient  stratification  with  respect  to  tumour  oxygenation 
status  could  be  clinically  important  (Hockel  et  al  1996,  Fyles  et  al  1998,  Welch  et  al  2003).  It 
has  been  hoped  that  modulation  of  tumour  oxygenation  could  be  applied  to  enhance  therapeutic 
efficacy.  An  attractive  intervention  is  breathing  hyperoxic  gas,  and,  indeed,  several  clinical 
trials  have  examined  the  efficacy  of  normobaric  or  hyperbaric  oxygen,  to  improve  therapeutic 
outcome,  but  often  with  marginal  success  (Overgaard  and  Horsman  1996).  It  has  been 
suggested  that  outcome  might  have  been  improved,  if  responsive  tumours  could  have  been 
identified  a  priori.  Accordingly,  accurate  evaluation  of  tumour  oxygenation  in  response  to 
interventions  at  various  stages  of  growth  should  provide  a  better  understanding  of  tumour 
response  to  therapy,  potentially  allowing  therapy  to  be  tailored  to  individual  characteristics. 

Given  the  importance  of  tumour  oxygenation,  many  techniques  have  been  developed  based 
on  microelectrodes,  optical  reflectance,  electron  paramagnetic  resonance  (EPR),  magnetic 
resonance  imaging  (MRI)  and  nuclear  medicine  approaches,  as  reviewed  previously  (Mason 
et  al  2002,  Zhao  et  al  2004).  While  each  approach  has  unique  strengths,  some  are  highly 
invasive.  Near-infrared  spectroscopy  (NIRS)  has  been  developed  in  recent  years  as  a 
promising  non-invasive  technique  to  quantify  the  concentration  of  tissue  chromophores,  such 
as  oxygenated  and  deoxygenated  haemoglobin,  water  and  lipid  (Sevick  et  al  1991,  Liu  et  al 
2000).  Due  to  the  deep  penetration  depth  and  biochemical  specificity  of  NIRS,  it  has  been 
widely  applied  for  quantitative  measurements  of  cerebral  oxygenation  (Delpy  and  Cope  1997, 
Yodh  and  Boas  2003)  and  blood  oxygenation  in  muscles  in  vivo  (Homma  etal  1996).  Recently, 
NIRS  has  been  also  used  to  monitor  tumour  vascular  oxygenation  with  respect  to  interventions 
(Hull  et  al  1999,  Liu  et  al  2000,  Gu  et  al  2003).  With  notable  exceptions  (Chance  1997, 
Dehghani  et  al  2004,  Shah  et  al  2004,  Pogue  et  al  2003),  NIRS  currently  lacks  spatial 
resolution,  and,  thus,  the  utility  of  global  measurements  requires  validation,  given  the  well- 
documented  heterogeneity  of  tumour  oxygenation.  In  this  regard,  Conover  et  al  (2000) 
compared  the  spatially  averaged  measurement  of  tumour  oxygen  saturation  (SOi)  using  NIRS 
with  the  local  SO2  in  individual  blood  vessels  measured  by  cryospectrophotometry.  The 
sensitivity  and  specificity  analysis  suggests  that  NIRS  may  identify  clinically  relevant  hypoxia, 
even  when  its  spatial  extent  is  below  the  resolution  limit  of  the  NIRS  technique.  We  have 
previously  investigated  correlates  between  p02  assessed  by  electrodes  (Mason  et  al  2003, 
Kim  et  al  2003b)  or  fibre-optic  probes  and  NIRS  (Gu  et  al  2003).  On  occasion,  there  was  a 
good  correlation  between  global  vascular  oxygenation  and  local  p02  at  individual  locations, 
but  often,  disparate  behaviour  was  observed.  Sequential  MRI  and  NIRS  suggested  a  better 
relationship  based  on  average  p02  from  multiple  locations  (Kim  et  al  2003b).  We  have  now 
implemented  simultaneous  NIRS  and  19F  MRI  to  examine  the  relationships  further.  Since 
NIRS  is  entirely  non-invasive  it  would  provide  an  attractive  surrogate  for  monitoring  tumour 
oxygenation,  and,  hence,  we  seek  correlations  with  absolute  p02  measurements  observed 
simultaneously  by  MRI. 

2.  Materials  and  methods 

Investigations  were  approved  by  the  Institutional  Animal  Care  and  Use  Committee. 

2.1.  Animal  preparation  and  experimental  set-up 

Mammary  adenocarcinomas  13762NF  (cells  originally  obtained  from  the  Division  of  Cancer 
Therapeutics,  NCI)  were  implanted  in  skin  pedicles  (Hahn  et  al  1993)  on  the  foreback  of  ten 
adult  female  Fisher  344  rats  ( — 150  g).  When  the  tumours  reached  ~1  cm  in  diameter,  the  rats 
were  anaesthetized  with  ketamine  hydrochloride  i.p.  (100  mg  kg-1  body  weight,  Aveco,  Fort 
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Figure  1.  Schematic  of  experimental  set-up.  Z  is  along  the  bore  of  the  magnet,  and  X  along 
the  axis  of  the  RF  coil.  PMT  represents  a  photomultiplier  tube.  I/Q  represents  an  In-phase  and 
quadrature  (I/Q)  demodulator  chip  and  ADC  represents  an  analogue-to-digital  converter. 


Dodge,  IA)  and  were  maintained  under  general  gaseous  anaesthesia  (air  and  1%  isoflurane; 
Baxter  International  Inc,  Deerfield,  IL).  Tumour  hair  was  trimmed  to  give  good  optical  contact 
for  NIR  light  transmission.  Hexafluorobenzene  (HFB,  50  /xl ,  99.9%,  Lancaster  Co,  Pelham, 
NH)  was  administered  along  two  or  three  tracks  in  central  and  peripheral  regions  of  the  tumours 
in  a  single  plane  (transverse  to  the  rat’s  tumour,  and  in  the  region  of  NIR  photon  pathway) 
using  a  Hamilton  syringe  with  a  32  G  needle.  The  needle  was  inserted  manually  to  penetrate 
across  the  tumour  and  was  withdrawn  ~1  mm  to  reduce  pressure  and  3-4  //I  of  HFB  were 
deposited.  The  needle  was  then  repeatedly  withdrawn  1-2  mm  and  further  HFB  deposited  at 
each  point,  as  described  in  detail  previously  (Zhao  et  al  2004). 

The  tumour  was  placed  inside  a  size-matched  Helmholtz  coil,  specially  designed  for 
the  simultaneous  MRI-NIRS  study.  The  tumour  was  inserted  between  the  two  loops  of  the 
Helmholtz  coil  and  two  NIRS  probes  were  introduced  through  the  ends  of  the  coil  along  the 
coil  axis  (figure  1).  The  probes  were  positioned  to  be  in  the  same  plane  as  the  HFB  injection. 
The  rats  were  placed  in  the  magnet  on  their  side,  and  body  temperature  was  maintained 
using  a  warm  water  blanket.  A  total  of  ten  rats  were  used  in  the  study:  seven  rats  were 
subjected  to  respiratory  challenge  in  the  sequence  of  air-oxygen-air-carbogen-air,  one  rat 
breathed  air-carbogen-air-oxygen-air,  one  rat  breathed  air-carbogen-air,  and  one  breathed 
air-oxygen-air. 

2.2.  NIRS  for  measuring  A[Hb02] 

A  homodyne  frequency-domain  system  (NIM,  Philadelphia,  PA)  was  used  to  monitor  the  global 
change  of  deoxy-  and  oxyhaemoglobin  concentration  (AlHbCL])  in  the  tumour,  as  described 
previously  (Yang  et  al  1997),  though  with  minor  modifications  to  ensure  MR  compatibility 
(figure  1).  Briefly,  the  laser  light,  which  was  emitted  from  two  laser  diodes  (at  758  nm  and 
785  nm),  was  amplitude-modulated  at  140  MHz  and  time  gated  on  and  off  sequentially.  The 
two  time-shared  laser  beams  illuminated  the  tumour  surface  alternately  through  a  light-delivery 
fibre  bundle  (source  fibre)  with  a  7  m  length  and  a  3  mm  bundle  diameter.  The  long  fibre 
bundle  ensured  the  separation  of  the  NIRS  hardware  from  the  magnet.  After  being  absorbed 
and  scattered  in  the  tumour  tissue,  the  transmitted  light  was  collected  on  the  opposite  side  of  the 
tumour  by  another  fibre  bundle  (detector  fibre)  of  same  length  and  diameter  as  those  of  source 
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fibre.  The  collected  optical  signal  was  detected  and  amplified  by  a  photomultiplier  tube  (PMT). 
The  probe  sheaths  at  the  fibre  tips  were  made  of  nylon  for  MR  compatibility.  An  In-phase  and 
quadrature  (I/Q)  demodulator  chip  was  used  to  demodulate  the  amplitude-modulated  signal 
from  the  PMT. 

In  principle,  since  the  I/Q  system  could  give  both  phase  and  amplitude  values,  it  should 
be  possible  to  calculate  absolute  HbC>2,  Hb  and  SO2  (Yang  et  al  1997).  However,  given  the 
tumour’s  small  size  and  large  spatial  heterogeneity,  we  are  currently  unable  to  achieve  reliable 
absolute  quantification  using  the  conventional  algorithm  (Fishkin  and  Gratton  1993)  based  on 
the  diffusion  approximation.  Instead,  based  on  modified  Beer-Lambert’s  law,  we  simply  use 
the  amplitude  of  the  light  transmitted  through  the  tumour  to  calculate  concentration  changes 
in  HbCT,  Hb  and  Hbt  of  the  tumour  caused  by  respiratory  intervention.  Namely,  changes  of 
oxy-  and  deoxyhaemoglobin  concentration,  A  [HbCT]  and  A[Hb],  respectively,  can  be  derived 
from  the  measured  amplitudes  at  the  two  wavelengths  (758  nm  and  785  nm)  (Kim  et  al  2003b) 
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where  Aj,  is  the  baseline  amplitude.  A,  is  the  transient  amplitude  during  measurement  and 
d  is  the  direct  source-detector  separation.  DPF  (differential  path-length  factor)  is  a  tissue- 
dependent  parameter  and  is  defined  as  the  ratio  between  the  optical  path  length  and  the  physical 
separation  between  the  source  and  detector.  The  units  for  AfHbCy,  A[Hb]  and  A[Hb]lolai  are 
mM.  DPF  is  a  variable  depending  on  tissue  type  and  wavelengths  and  it  is  currently  difficult  to 
quantify  DPF  for  tumours.  Since  this  study  focuses  on  dynamic  changes  of  [HbCT],  A[HbCT] 
values  may  be  scaled  by  a  factor  of  DPF  (in  units  of  mM/DPF),  to  obtain  characteristic  features 
of  tumour  oxygen  dynamics  (Liu  et  al  2000). 


2.3.  Mathematical  model  for  blood  oxygenation  dynamics  of  tumours 


Based  on  our  previous  experimental  study,  we  derived  a  simple  mathematical  model  (Liu 
et  al  2000)  to  examine  tumour  vascular  dynamics  during  oxygen  intervention,  by  analogy  to 
the  method  used  to  quantify  regional  cerebral  blood  flow  (rCBF)  with  diffusible  radiotracers, 
originally  developed  by  Kety  (1951).  Accordingly,  changes  of  oxyhaemoglobin  concentration 
in  tumour  vasculature  induced  by  hyperoxic  gas  intervention,  A[HbCT],  can  be  expressed  as 


A[Hb02(r)rsculamre 


yH0 
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where  A  [HbCT]  corresponds  to  the  changes  in  oxyhaemoglobin  concentration  from  tumour 
vasculature  measured  by  NIRS,  H0  is  the  arterial  oxygenation  input  for  AHbCTaltery  after 
time  0,  /  is  the  blood  perfusion  rate  and  y  is  the  vascular  coefficient  of  the  tumour,  defined 
as  the  ratio  of  [HbCT]  changes  in  the  vascular  bed  to  that  in  the  veins.  Specifically,  y  = 
A[Hb02]vasculature/A[Hb02]vem,  r  is  the  time  constant  ( =y/f )  and  A  =  yH0. 

Given  that  solid  tumours  often  develop  hypoxic  regions,  which  are  poorly  perfused  (Jiang 
et  al  2004,  Mason  et  al  1994,  Mazurchuk  et  al  1999,  Song  et  al  2002),  we  hypothesized 
that  the  observed,  bi-phasic  feature  of  A[Hb02]  came  from  two  different  perfusion  regions 
(well-  and  poorly  perfused  regions).  Therefore,  it  is  reasonable  to  define  two  different  blood 
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perfusion  rates  (/),  f2)  with  two  different  vascular  coefficients  (y i,  y2)  in  the  mathematical 
model.  Consequently,  equation  (3)  has  been  modified  to  a  summation  of  two  exponential 
expressions,  representing  two  perfusion  regions,  as 
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where  f\  and  y  \  are  the  blood  perfusion  rate  and  the  vasculature  coefficient,  respectively,  in 
the  well-perfused  region,  and  f2  and  y2  represent  the  same  respective  meanings  in  the  poorly 
perfused  region.  Also,  it  follows  that  Ai  =  Y\Hq,  A2  =  y2Ho,  x\  =  yi/fi  and  x2  =  y2/  f2. 
Since  A\,  A2,  r  1  and  t2  can  be  determined  by  curve-fitting  equation  (4)  to  the  dynamic  NIRS 
measurements,  we  can  obtain  the  ratios  for  two  vascular  coefficients  and  two  blood  perfusion 
rates,  as  follows  (Liu  et  al  2000): 


n  _  An  fi  _  A\/r\ 

y2  A2  f2  A2/x2 


These  two  ratios  provide  insight  into  tumour  vascular  structures  and  blood  perfusion  rates. 
Specifically,  y\/y2  may  be  associated  with  the  vascular  volume  fraction  between  the  two 
regions,  and  /1//2  reflects  the  perfusion  ratio  between  the  two  regions.  For  example,  if  /1//2 
is  substantially  greater  than  1,  the  two  compartments  have  significantly  different  perfusion 
rates. 

This  mathematical  model  is  consistent  with  experimental  animal  studies  (Liu  et  al  2000, 
Gu  et  al  2003,  Kim  et  al  2003b).  Moreover,  it  is  supported  by  computer  simulations  (Kim  and 
Liu  2005)  and  a  phantom  study  (Kim  and  Liu  2004)  each  of  which  shows  that  the  bi-phasic 
or  bi-exponential  feature  can  be  present  if  both  a  slow  and  a  fast  perfusion  of  flow  co-exist 
within  the  interrogated  volume. 


2.4.  FREDOM  for  measuring  p02 

MRI  was  performed  using  a  Varian  Inova  4.7  T  horizontal  bore  system  equipped  with  actively 
shielded  gradients.  Shimming  was  performed  on  the  tumour  tissue  water  signal  to  reduce 
the  line- width  to  less  than  100  FIz.  'H  MRI  (200.1  MHz)  T1 -weighted  reference  images 
were  acquired  with  TR/TE  of  150/10  ms  and  40  x  40  mm2  field  of  view.  Following  ]H 
MRI,  corresponding  19F  MR  images  (188.3  MHz)  with  a  matrix  size  of  32  x  32  (1.25  mm 
per  pixel)  were  obtained  to  show  the  distribution  of  HFB  in  the  tumour.  A  single  10  mm 
thick  slice  was  chosen  to  include  all  the  injected  HFB,  as  the  plane  of  injection  was  not 
always  perfectly  aligned  with  the  imaging  plane.  The  actual  volume  interrogated  is  defined 
by  the  distribution  of  the  reporter  molecule  in  the  third  dimension  rather  than  the  imaging 
gradients  (Zhao  et  al  2004).  The  FREDOM  (Fluorocarbon  A’elaxometry  by  Fcho-planar 
imaging  for  Dynamic  Oxygen  Mapping)  approach  was  used  to  measure  pCL,  as  described  in 
detail  previously  (Hunjan  et  al  2001).  The  spin  lattice  relaxation  rate  (R1  =  1/T1 )  of  HFB  is 
highly  sensitive  to  changes  in  p02,  but  not  sensitive  to  temperature  variations:  a  deviation  of 
1  °C  in  temperature  will  introduce  a  deviation  of  only  0.13  Torr  in  the  p02  estimate,  when  p02 
is  about  5  Torr  (Zhao  et  al  2004).  T1  maps  were  computed  on  a  voxel-by-voxel  basis  using 
nonlinear  least-squares  data  fitting  by  the  Gauss-Newton  method.  We  applied  a  threshold  to 
the  raw  T1  data  in  order  to  remove  random  noise,  i.e.,  voxels  with  T1  error  >3.6  s  or  T1 
error/Tl  >50%  were  disregarded.  Maps  of  p02  values  with  1.25  mm  pixel  size  were  obtained 
from  the  T1  maps  using  the  equation,  p02  =  (1/T1  —  0.0835)/0.001  876  (Hunjan  et  al  2001). 
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Figure  2.  Temporal  profiles  of  A[[Ib(T|  (curve  only)  and  mean  p(T  (curve  with  ♦)  in  response  to 
respiratory  challenge.  Both  A[HbOi]  and  pC>2  are  group-averaged  data  from  seven  13762NF 
rat  breast  tumours  with  an  intervention  sequence  of  air-oxygen-air-carbogen-air,  measured 
simultaneously  by  NIRS  and  FREDOM.  The  standard  deviations  of  the  volume-averaged  p02 
from  seven  tumours  are  indicated  for  baseline,  but  became  very  large  in  response  to  interventions 
and  are  omitted.  AlHbOi]  has  units  of  mM/DPF. 


Each  pCE  map  was  acquired  in  min.  Three  baseline  pC>2  data  sets  were  acquired  over 
24  min  for  all  tumours,  while  the  rats  breathed  air,  after  which  the  rats  were  repeatedly  exposed 
to  oxygen  or  carbogen  (95%  O2  and  5%  CO2)  interventions.  Five  pCB  maps  were  obtained 
during  each  subsequent  gas  switch  period.  Typically,  for  a  five  gas-intervention  sequence 
(e.g.,  air-oxygen-air-carbogen-air),  a  total  of  23  pCB  maps  were  obtained  over  a  period  of 
3  h.  Due  to  thresholding,  some  voxels  did  not  appear  in  all  pCT  maps  and  these  were  discarded. 
For  temporal  analysis,  voxels  were  selected  as  only  those  which  provided  consistently  reliable 
data  (i.e.,  satisfied  all  the  thresholding  criteria  specified  above)  for  all  23  measurements  over 
the  time  course  with  a  range  of  5 — 44  acceptable  voxels  for  the  ten  tumours.  The  slope  of 
dynamic  pC>2  changes  (rate)  was  defined  as  d(pC>2)/dr  and  dfpCBO/df  in  response  to  increasing 
or  decreasing  inhaled  FO2  (Fraction  of  O2),  respectively. 

2.5.  Statistical  analysis 

Finear  regression  analysis  was  used  to  calculate  the  correlation  between  the  NIRS-derived 
tumour  haemodynamic  parameters  (i.e.,  AfHbCB],  A\/t\,  A2/T2,  and  the  FREDOM- 

determined  tumour  parameters  (i.e.,  pOy,  d(pC>2)/df,  AipOfl/dt).  Data  are  presented  as  mean 
±  standard  deviation  (SD)  and  paired  Student-t  tests  compared  the  effects  of  oxygen  and 
carbogen  on  A[Hb02]  and  pCK 


3.  Results 

3.1.  Dynamic  response  of  N[HbO  2]  measured  by  NIRS 

Figure  2  shows  a  group-averaged,  temporal  profile  of  A[Hb02]  from  seven  13762NF 
rat  breast  tumours  with  air-oxygen-air-carbogen-air  intervention,  displaying  apparent 
biphasic  responses  to  both  interventions.  Both  single  (equation  (3))  and  double-exponential 
(equation  (4))  curve  fitting  were  tested  for  the  carbogen  intervention  in  one  representative 
tumour  (1.6  cm3),  as  in  figure  3(A).  The  maximal  AfHbCT]  achieved  with  oxygen  challenge 
was  compared  with  that  of  carbogen,  and  revealed  no  significant  difference  between  oxygen 
and  carbogen  interventions  (p  >  0.3);  indeed,  there  was  a  strong  correlation  between  the 
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Figure  3.  (A)  Dynamic  responses  of  A[HbC>2]  (rnM/DPF)  to  gas  intervention  (carbogen).  Single 
exponential  curve  fitting  yielded  A[Hb(>>]  =  0.22(1  —  exp[— (t  —  106.8)/2.5]}  (R2  =  0.64),  and 
double  exponential  fitting  resulted  in  A 1 1  IbOi |  =  0.15(1  —  exp[— (r  —  106.8)/0.58]}+  0.13(1  — 
exp[— (f  —  106.8)/23.2])  (R1  =  0.88).  (B)  Relationship  of  maximum  A[HbOi]  (mM/DPF)  in 
breast  tumours  in  response  to  switching  from  air  to  oxygen  and  to  carbogen  (R2  >  0.75). 

maximal  A[HbOi]  values  with  these  two  interventions  ( R 2  >  0.75,  figure  3(B)),  consistent 
with  our  previous  observations  (Gu  et  al  2003).  No  correlation  was  found  between  the 
perfusion  rate  ratio  (/1//2)  and  tumour  size  ( R 1  =  0.16).  Vascular  oxygen  dynamics  in 
response  to  interventions  are  provided  for  individual  tumours  in  table  2. 

3.2.  pC>2  measurements  by  FREDOM 

Overlay  of  19F  and  'FI  MR  images  demonstrated  that  HFB  was  distributed  in  both  central 
and  peripheral  regions  of  tumours  (not  shown).  Individual  pOy  values  taken  from  each  voxel 
ranged  from  hypoxic  (<  1  Torr)  to  35  Torr  under  baseline  conditions  (figure  4).  Mean  baseline 
PO2  values,  which  are  averaged  over  all  voxels  in  a  given  tumour  in  the  first  three  p02 
maps,  ranged  from  hypoxic  (<5  Torr)  to  27  Torr  with  a  hypoxic  fraction  (HF5;  fractional 
voxels  that  are  less  than  5  Torr)  ranging  from  0  to  100%  (mean  36%)  and  summarized 
in  table  1.  A  strong  correlation  was  found  between  mean  baseline  pC>2  and  HF5  (R1  >  0.85, 
figure  5).  Administration  of  oxygen  or  carbogen  produced  significant  increases  in  tumour  pC>2, 
as  shown  in  the  pCF  maps,  graphs  and  table.  The  tumour  pCF  values  were  averaged  over  the 
entire  slice  and  responses  to  respiratory  challenge  for  the  group  of  seven  tumours,  measured 
simultaneously  by  FREDOM  and  NIRS,  are  shown  in  figure  2.  Baseline  measurements 
(breathing  air)  were  generally  stable,  and  altering  the  inhaled  gas  to  oxygen  or  carbogen 


Table  1.  Oxygen  tension  (PO2)  in  ten  rat  mammary  13762NF  adenocarcinomas.  p02  in  baseline  is  the  average  value  of  all  voxels  in  the  first  three  maps,  while  p02  for  oxygen  or 
carbogen  is  the  average  value  of  all  voxels  in  the  final  three  p02  maps  during  exposure  to  the  hyperoxic  gas.  Hypoxic  fraction  (HF)  is  the  percentage  of  voxels  with  p02  values  less  than 
5  Torr  (HF5)  or  10  Torr  (HF10)  and  HF5  and  HF10  are  the  average  values  of  final  three  maps  in  baseline,  oxygen  or  carbogen  respectively.  The  mean  p02  increased  significantly  with 
both  oxygen  and  carbogen  and  the  value  was  significantly  higher  with  carbogen.  HF5  was  significantly  reduced  with  oxygen  or  carbogen. 


Tumour 

Baseline  (air) 

Oxygen  challenge 

Carbogen  challenge 

No 

Volume  (cm3) 

Voxels3 

pCh  (Torr) 
Mean  ±  SD 

HF5  (%) 

HF10  (%) 

pCF  (Torr) 
Mean  ±  SD 

HF5  (%) 

HFio  (%) 

pOj  (Torr) 
Mean  ±  SD 

HF5  (%) 

HFio  (%) 

1 

1.6 

26 

12  ±3 

24  ±  14 

49  ±  15 

56  ±  10 

22  ±  2 

29  ±8 

66  ±3 

0 

0 

2 

0.3 

10 

12±  2 

6±  12 

53  ±  15 

62  ±  15 

0 

0 

79  ±6 

0 

0 

3 

0.5 

5 

4  ±  4 

60  ±20 

80  ±0 

14  ±  1 

33  ±  12 

53  ±  12 

3  ±  1 

73  ±  12 

73  ±  12 

4 

1.3 

14 

<1  Torrb 

100  ±0 

100  ±  0 

14  ±7 

19  ±  1 1 

48  ±  30 

61  ±  21 

0 

0 

5 

0.9 

24 

24  ±4 

0  ±  0 

8  ±  4 

47  ±  10 

0 

0 

71  ±  11 

0 

0 

6 

1.2 

8 

4  ±  1 

58  ±7 

63  ±0 

43  ±  12 

50  ±  22 

63  ±0 

73  ±  10 

0 

0 

7 

0.8 

24 

27  ±2 

0  ±  0 

3  ±  2 

122  ±  21 

0 

0 

142  ±  26 

0 

0 

8 

0.6 

8 

9  ±  6 

50  ±22 

63  ±22 

103  ±  26 

0 

0 

118  ±37 

0 

0 

9 

0.7 

44 

25  ±  1 

10  ±5 

26  ±3 

- 

- 

- 

35  ±3 

0 

3  ±  3 

10 

0.7 

8 

7  ±  6 

54  ±  19 

75  ±  22 

52  ±  18 

0 

0 

- 

- 

- 

Mean  ±  SD  0.9  ±  0.4 

17  ±  12 

12  ±  10 

36  ±33 

52  ±31 

57  ±  36° 

14  ±  ld 

21  ±27 

72  ±  41b 

8  ±  24b 

8  ±24 

a  Number  of  voxels  that  provide  acceptable  data  for  all  23  pC>2  maps. 

bp  <  0.01. 

cp<  0.005. 

Ap  <0.05. 
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Figure  4.  p02  maps  (1.25  mm  resolution)  obtained  using  FREDOM  overlaid  on  the  1 1  anatomic 
image  of  a  tumour  (FOV  =  4  cm  x  4  cm).  (Left)  rat  breathing  air,  (centre)  breathing  100%  O2 
(fifth  map  after  switching  from  air),  (right)  breathing  carbogen  (95%  O2  +  5%  CO2)  (fifth  map 
after  switching  from  air).  The  pC>2  maps  show  distinct  heterogeneity. 

(This  figure  is  in  colour  only  in  the  electronic  version) 


Baseline  HF5  (%) 

Figure  5.  Correlation  between  baseline  p02  and  hypoxic  fraction  (HF5)  measured  using  FREDOM 
(R2  >  0.85). 

induced  rapid  and  significant  changes  in  both  p02  and  A[Hb02]  (p  <  0.001)  Upon  return 
to  air  (baseline),  A[Hb02]  dropped  quickly  and  significantly  within  16  min,  and  then  more 
slowly,  for  the  next  24  min,  whereas  the  p02  decrease  was  more  gradual.  Altering  the  inhaled 
gas  to  carbogen  also  produced  a  rapid  increase  in  both  p02  and  A[Hb02].  Upon  return  to 
air  breathing  from  carbogen,  both  A[Hb02]  and  p02  showed  a  similar  trend  to  that  following 
oxygen.  As  expected,  all  ten  tumours  showed  a  significant  increase  in  p02  (average  over  entire 
slice),  and  decrease  in  hypoxic  fraction  (HF)  in  response  to  oxygen  or  carbogen  inhalation. 
The  magnitude  of  response  to  either  hyperoxic  gas  was  correlated  ( R 2  >  0.79),  as  was  the 
maximum  volume-averaged  p02  achieved  with  either  gas  ( R 2  >0.83).  The  rate  of  increase  with 
oxygen  challenge,  d(p02)/df,  was  significantly  faster  than  the  return  to  baseline,  c/(p02')/df, 
for  oxygen  intervention  (p  <  0.02),  but  no  difference  was  observed  for  carbogen  ( p  >  0.1, 
figure  6).  The  mean  p02  values  of  individual  tumours  averaged  over  the  final  three  p02  maps 
during  exposure  to  carbogen  breathing  were  significantly  higher  than  oxygen  ( p  <  0.01),  and 
the  tumour  hypoxic  fraction  was  generally  eliminated  during  carbogen  breathing  (n  =  7  of  9 
tumours,  table  1). 

3.3.  The  relationship  between  p02  and  A[HbC>2]  with  respect  to  hyperoxic  gas 

Taken  as  a  group  of  ten  tumours,  there  was  no  apparent  relationship  between  the  magnitude 
of  the  change  in  tumour  vascular  oxygenation  (A[Hb02])  and  change  in  p02  (R2  <  0.1). 
However,  if  tumours  were  divided  into  two  sub-populations,  then  two  separate  correlations 
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Figure  6.  Mean  d(pC>2)/dr  (open)  and  d(pC>2 ')/dt  (shaded)  ±  SD  for  eight  tumours  with  both 
interventions,  when  gas  was  switched  from  air  to  the  hyperoxic  gas  and  back  to  air,  respectively. 
PO2  is  the  mean  value  of  all  acceptable  voxels  appearing  in  the  five  maps  during  the  oxygen  or 
carbogen  intervention.  d(pC>2)/dr  is  the  slope  of  regression  line  of  five  pC>2  readings  versus  time 
when  switching  from  air  to  oxygen  or  carbogen,  and  d(p02')/cfi  is  the  slope  of  regression  line  for 
five  p02  readings  during  the  switch  back  to  air.  The  rates  showed  a  significant  difference  with 
oxygen  ( p  <  0.02),  but  not  with  carbogen  ( p  >  0.1). 


0  0.1  0.2  0.3  0.4  0.5  0.6 

A[Hb02] 

(A) 


(B) 


Figure  7.  (A)  Correlation  between  maximum  A[HbC>2]  and  change  in  pC>2  with  respect  to 
hyperoxic  gas  intervention  for  two  groups  of  tumours  (•)  (group  1:  no  1,  2,  6,  7,  8,  9,  10;  R2  > 
0.51)  and  (A)  (group  2:  no  3, 4,  5;  R2  >  0.82).  The  unit  for  A[HbC>2]  is  mM/DPF.  (B)  Correlation 
between  mean  pC>2  achieved  with  hyperoxic  gas  breathing  and  perfusion  rate  ratio  (/1//2)  for 
tumours  with  biphasic  response  to  intervention  ( R 2  >  0.7).  pC>2  is  the  mean  value  for  the  final 
three  pC>2  maps  under  hyperoxic  intervention,  selected  from  all  voxels  appearing  in  the  final  three 
PO2  maps  during  oxygen  (♦)  or  carbogen  (□)  intervention. 


were  found  each  with  similar  slope  (figure  7(A)).  There  was  also  a  correlation  (. R 2  >  0.7) 
between  the  perfusion  rate  ratio,  /1//2,  derived  from  fitting  the  A[Hb02]  curve  and  the  mean 
pC>2  values  achieved  with  hyperoxic  gas  intervention  (figure  7(B)).  Assessment  of  /1//2  is 
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Figure  8.  d(pC>2)/dt  versus  A  \/z\  (determined  from  AfHbCh]  for  tumours  with  biphasic  response 
to  interventions)  showing  a  positive  correlation  ( R 2  >  0.5).  pC>2  values  are  the  mean  value  of  all 
well-oxygenated  voxels  (i.e.,  the  maximum  pC>2  >10  Torr  under  oxygen  or  carbogen  intervention) 
appearing  in  the  five  pC>2  maps  during  oxygen  (♦)  or  carbogen  (□)  intervention. 


Table  2.  Summary  of  vascular  oxygen  dynamics  for  the  experimental  tumours.  Two  amplitudes 
(Ai,  A2)  and  two  time  constants  (ti,  T2)  were  determined  by  curve-fitting  the  dynamic  NIRS 
measurements  using  a  double-exponential  expression.  Nine  out  of  ten  tumours  were  observed 
to  have  double-exponential  features  with  either  oxygen  or  carbogen  intervention.  A\  is 
significantly  smaller  than  A  2  with  oxygen  ( p  <  0.01),  but  no  significant  differences  in  carbogen 
(p  >  0.19). 


Tumour 

Double  exponential  fitting 

No 

Intervention 

A!  (AU) 

A2  (AU) 

z\  (min) 

12  (min) 

R 

la 

Oxygen 

0.11  ±  0.002 

0.12  ±  0.002 

0.49  ±  0.02 

6.6  ±  0.2 

0.94 

Carbogen 

0.15  ±0.001 

0.13  ±0.002 

0.58  ±0.01 

23.2  ±  1.1 

0.94 

2a 

Oxygen 

0.03  ±  0.002 

0.11  ±0.005 

0.46  ±0.18 

25.6  ±  2.9 

0.91 

Carbogen 

- 

- 

- 

- 

- 

3a 

Oxygen 

0.15  ±0.001 

0.38  ±  0.003 

0.31  ±0.02 

20.1  ±0.39 

0.98 

Carbogen 

0.23  ±  0.005 

0.2  ±  0.004 

1.2  ±0.04 

10.8  ±  0.45 

0.95 

4a 

Oxygen 

- 

- 

- 

- 

- 

Carbogen 

0.16  ±  0.002 

2.3  ±  0.27 

0.29  ±  0.02 

155.8  ±21 

0.97 

5a 

Oxygen 

0.12  ±  0.002 

0.23  ±  0.002 

0.55  ±  0.03 

11.6  ±0.27 

0.94 

Carbogen 

- 

- 

- 

- 

- 

6a 

Oxygen 

0.06  ±0.001 

0.26  ±  0.005 

0.13  ±0.01 

31.2  ±  1.1 

0.96 

Carbogen 

0.03  ±  0.001 

0.08  ±0.001 

0.04  ±  0.02 

13.9  ±0.6 

0.78 

7a 

Oxygen 

- 

- 

- 

- 

- 

Carbogen 

- 

- 

- 

- 

- 

8b 

Oxygen 

0.11  ±  0.002 

0.2  ±  0.001 

0.24  ±  0.02 

7.1  ±  0.1 

0.96 

Carbogen 

0.2  ±  0.002 

0.12  ±  0.001 

0.49  ±  0.01 

9.1  ±0.2 

0.95 

9C 

Carbogen 

0.01  ±  0.0004 

0.02  ±  0.0004 

0.14  ±0.02 

3.9  ±0.12 

0.64 

10d 

Oxygen 

0.05  ±  0.002 

0.03  ±  0.002 

0.46  ±  0.05 

5.6  ±0.56 

0.68 

Mean  ±  SD 

0.11  ±  0.068 

0.32  ±  0.59 

0.42  ±  0.3 

24.9  ±  40.2 

- 

a  Air  — >  O2  — >■  air  — >  carbogen  — ►  air. 
b  Air  — ►  carbogen  — ►  air  — »  O2  — *  air. 
c  Air  — ►  carbogen  — ►  air. 
d  Air  — >  oxygen  — >  air. 


56 


M  Xia  el  al 


predicated  on  biphasic  behaviour  with  respect  to  interventions,  which  was  observed  in  most 
cases  (13  of  16  measurements).  There  was  also  apositive  correlation  between  Ai/tj  (the  fast 
component  of  biphasic  AfHbCT])  and  the  d(pC>2)/df  of  well-oxygenated  voxels  (i.e.,  those  with 
pCK  values  >10  Torr  under  oxygen  or  carbogen  intervention)  (R1  >  0.5,  figure  8).  However, 
no  correlation  was  found  between  d(pC>2)/dr  and  A2IT2  (the  slow  component). 

4.  Discussion 

Integration  of  diverse  imaging  techniques  can  be  technically  challenging,  since  each  modality 
has  specific  technical  constraints  and  requirements.  Here,  the  NIR  system  was  modified  to 
use  longer  optical  fibres  and  any  metal  components  were  eliminated  from  the  fibre  tips.  Due 
to  the  spatial  restrictions  within  the  bore  of  the  magnet  a  Helmholtz  coil  was  built  specifically 
providing  access  to  both  the  tumour  and  fibres.  The  fibres  required  sufficient  flexibility  to 
allow  them  to  be  bent  through  requisite  angles  within  the  confines  of  the  magnet  bore.  To  date 
there  have  been  few  reports  of  simultaneous  data  acquisition  by  NIR  and  MRI,  e.g.,  studies  of 
phantoms  (Pogue  et  al  2003),  human  brain  (Toronov  et  al  2001,  Chen  et  al  2003)  and  breast 
(Ntziachristos  et  al  2000,  Gulsen  et  al  2002,  Ntziachristos  et  al  2002,  Gu  et  al  2004). 

In  the  present  study,  global  average  A[Hb02]  was  measured  by  NIRS,  and  p02  maps 
were  obtained  simultaneously  by  19F  MRI.  We  used  transmission  mode  NIRS  in  order  to 
interrogate  deep  tumour  tissue.  While  NIRS  is  a  global  measurement,  the  region  sampled  by 
it  is  predominantly  a  banana-shaped  region  joining  the  locations  of  the  source  and  detector 
(Arridge  1995,  Arridge  and  Schweiger  1995),  and  we  recently  showed  that  typically  15  to 
30%  of  the  vascular  volumes  of  rat  tumours  are  interrogated  by  NIR  using  this  configuration 
(Gu  et  al  2005).  We  positioned  the  NIRS  probes  (source  and  detector)  to  be  in  the  same 
plane  as  the  HFB  injection  to  ensure  maximum  overlap  between  the  regions  sampled  by  the 
two  techniques.  Utilizing  our  previously  developed  mathematical  model  (Liu  et  al  2000), 
multiple  haemodynamic  parameters  were  derived  for  A[Hb02]  A2IX2  and  /1//2)  to 

be  compared  with  p(K  Our  results  demonstrate  that  oxygenation  parameters  measured  from 
both  techniques  show  significant  and  consistent  elevation  in  tumour  oxygenation  during  the 
hyperoxic  gas  interventions.  As  reported  previously,  the  magnitude  of  the  vascular  response 
was  similar  with  both  hyperoxic  gases  (Gu  et  al  2003).  As  expected,  A[Hb02]  increased 
much  faster  than  p02  in  all  ten  tumours,  indicating  that  change  in  tumour  vascular  oxygenation 
precedes  tumour  tissue  oxygenation.  This  observation  is  consistent  with  our  previous  studies 
in  this  tumour  type  measured  simultaneously  by  NIRS  and  fibre-optic  probes  (Gu  et  al  2003), 
as  well  in  the  Dunning  prostate  R3327-AT1  tumours  measured  sequentially  by  NIRS  and  19F 
MRI  (Kim  et  al  2003b). 

We  have  previously  demonstrated  the  application  of  FREDOM  to  monitor  tumour  oxygen 
dynamics  in  diverse  rat  prostate  tumours  (Zhao  et  al  2001,  2002),  human  tumour  xenografts 
(Mason  et  al  2002)  and  a  few  breast  tumours  (Song  et  al  2002).  Here,  the  mean  baseline  was 
PO2  =  12  ±  10  Torr  for  the  ten  tumours,  which  is  lower  than  reported  previously  (Song  et  al 
2002),  but  entirely  consistent  with  the  newer  anaesthetic  protocol  (air  or  21%  oxygen,  as 
opposed  to  33%  O2  previously).  There  was  a  strong  correlation  between  baseline  p02  and 
hypoxic  fraction  (figure  5),  as  we  have  previously  found  using  Dunning  prostate  R3327-HI 
tumours  (Zhao  et  al  2001).  The  pOa  achieved  with  carbogen  in  this  study  was  significantly 
higher  than  with  oxygen  and  carbogen  appeared  to  be  more  effective  at  eliminating  the  hypoxic 
fraction  (table  1).  However,  carbogen  was  generally  applied  second  in  our  experimental 
protocols,  and  it  is  highly  likely  that  the  initial  oxygen  primed  the  tumour.  Indeed,  while 
oxyhaemoglobin  generally  returned  to  baseline  during  the  air  breathing  episode  between 
hyperoxic  gases,  it  is  clear  that  pOi  remained  elevated  (figure  2).  Both  ApOi  and  the 
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maximum  pC>2  achieved  with  either  gas  were  closely  correlated.  Similar  behaviour  was 
reported  previously  based  on  measurements  using  fibre  optic  probes  (Gu  et  al  2003).  For 
the  group  of  tumours  in  this  study,  mean  baseline  pOy  did  not  provide  a  good  indication  of 
response  to  hyperoxic  gas  (ApCF  or  pCFmax)-  However,  considering  the  fate  of  individual 
voxels  in  separate  tumours  with  respect  to  intervention,  strong  linear  correlation  was  observed 
in  some  tumours  (tumour  no  3, 4,  6  and  9)  between  initial  mean  baseline  pCF  and  the  maximum 
pCF  at  the  same  location  (voxel)  in  response  to  carbogen  or  oxygen  breathing  consistent  with 
the  results  of  Song  et  al  (2002).  Other  tumours  showed  no  correlation. 

The  rate  of  p02  response  to  either  gas  was  similar  (figure  6),  but  decrease  upon  return  to 
air  was  significantly  faster  in  the  case  of  oxygen,  but  not  carbogen.  While  the  relationships 
between  Ap02  and  A|  HbOi  ]  are  not  obvious,  there  exists  a  significant  correlation  between 
f\  If 2  and  mean  p02  values  achieved  with  hyperoxic  gas  intervention.  Most  f\  / fi  values  are 
between  5  and  60  (figure  7(B)),  implying  that  the  blood  perfusion  rate  in  the  well-perfused 
region  is  much  higher  than  that  from  the  poorly  perfused  region.  Not  surprisingly  the  higher 
the  ratio  of  well  perfused  to  poorly  perfused  regions,  the  higher  the  mean  pCF  values  achieved 
by  the  hyperoxic  gas  interventions. 

Previous  studies  have  demonstrated  that  tumour  tissue  oxygenation  could  be  strongly 
affected  by  changes  in  tumour  blood  flow  measured  locally  by  laser  Doppler  flowmetry 
(Vaupel  et  al  1994).  It  was  also  reported  that  spontaneous  fluctuations  in  flow  and  peri¬ 
vascular  pCF  are  correlated  at  the  micro-regional  level  (Kimura  et  al  1996,  Braun  et  al  1999). 
In  our  study,  we  are  able  to  derive  the  perfusion  rate  ratios  between  the  well  perfused  and 
poorly  perfused  regions  from  the  haemoglobin  concentration  responses  to  the  intervention 
measured  by  NIRS.  We  found  that  there  was  a  significant  correlation  between  d(pC>2)/dr  and 
Ai/t  i  (/"],  proportional  to  Ai/rj),  but  not  d(p02)/d/  and  A2/ r2  (/o,  proportional  to  A2/r2), 
provided  that  the  pCF  readings  were  selected  from  well-oxygenated  or  responsive  voxels 
(figure  8).  This  linear  correlation  suggests  that  the  rate  of  change  in  pOy  is  closely  related 
to  the  perfusion  rate  in  the  well-perfused  region,  f\ .  In  other  words,  the  dynamic  changes  in 
pCF  of  those  regions  responsive  to  hyperoxic  gas  intervention  may  be  attributed  to  fast  tumour 
vascular  perfusion,  rather  than  to  the  slow  perfusion  in  tumour  vasculature.  We  believe  these 
results  provide  a  valuable  association  between  tumour  vascular  oxygenation  and  tumour  pCF 
determined  simultaneously  by  the  optical  and  NMR  measurements. 

A  goal  had  been  to  develop  a  low-cost,  simple,  fast  surrogate  measurement  of  pCF  based  on 
NIRS  of  the  oxygenation  status  of  endogenous  haemoglobin.  Both  figures  7(B)  and  8  suggest 
that  there  exist  linear  relationships  of  pCF  with  the  NIRS  measurable  parameters  f\  //>  and 
with  A\/x\.  However,  the  correlation  between  ApC>2  and  AlHbCB]  (figure  7(A))  seems  to 
be  more  complex  with  a  separation  of  the  tumours  into  two  groups.  It  may  be  noteworthy 
that  the  majority  (6  of  7)  of  the  tumours  associated  with  the  correlation  on  the  left-hand 
side  of  the  graph  had  a  high  initial  global  pCF  (i.e.,  mean  baseline  pCF  >5  Torr),  whereas 
2  out  of  3  of  the  tumours  on  the  right  had  low  initial  global  pCF.  A  possible  interpretation 
relates  to  the  shape  of  the  haemoglobin-oxygen  dissociation  curve  (the  Hill  curve).  A  given 
pCF  response  can  correspond  to  different  changes  in  haemoglobin  saturation  depending  on 
where  the  change  occurs  on  the  Hill  curve.  At  lower  initial  pCF  there  may  be  a  substantial 
ApCB  with  little  A  [HbOi].  By  contrast,  at  higher  initial  pCF  the  same  increase  in  pCF  could 
produce  a  greater  effect  on  the  saturation.  While  A  [HbOi]  is  proportional  to  the  change  of  SO2 
(with  the  assumption  of  little  change  in  total  haemoglobin  concentration  during  intervention), 
two  subpopulations  of  tumours  in  the  relationship  between  ApC>2  and  AfHbCB]  could  be  the 
result  of  different  initial  pCB.  Similarly,  two  subpopulations  of  tumours  were  observed  in  the 
relationship  between  initial  SO2  and  the  carbogen-induced  change  in  saturation  in  the  study  of 
Hull  et  al{  1 999).  Such  an  effect  also  confounds  the  direct  correlation  of  BOLD  MRI  response 
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to  changes  in  p02  (Baudelet  and  Gallez  2002).  Significantly,  preliminary  data  reported  by 
Gu  et  al  (2004)  for  simultaneous  BOLD  MRI  and  NIR  in  tumours  showed  a  strong  linear 
correlation  in  response  to  hyperoxic  gas  challenge,  as  also  reported  by  Chen  et  al  (2003)  in 
the  rat  brain. 

A  major  concern  is  tumour  heterogeneity,  as  recognized  throughout  the  literature  and 
shown  here  by  19F  MRI  (figure  4).  Indeed,  we  have  obtained  some  preliminary  data  using 
a  single  NIR  source  and  three  detectors  placed  on  various  regions  across  a  tumour  (Kim 
et  al  2003a),  showing  that  each  region  of  the  tumour  responded  differently  to  hyperoxic  gas, 
in  terms  of  the  extent  and  rate,  indicating  the  heterogeneity  of  tumour  vasculature.  Spatial 
discrimination  will  be  even  more  critical,  if  such  studies  are  transferred  to  human  breast 
cancer,  where  the  tumour  is  surrounded  by  normal  tissue  (Brooksby  et  al  2003).  Nevertheless, 
we  believe  this  haemodynamic  model  and  correlation  between  tumour  vascular  oxygenation 
and  p02  provides  valuable  insight  into  the  tumour  compartment  of  such  a  mixed  system  and 
explores  dynamic  signatures  of  breast  tumours,  which  could,  in  turn,  enhance/assist  human 
breast  cancer  diagnosis  and  prognosis. 

In  summary,  by  studying  tumour  vascular  oxygenation  concomitantly  with  changes  in 
tumour  oxygen  tension,  we  found  several  significant  correlations  between  rates  and  magnitudes 
of  vascular  and  tissue  responses.  This  study  also  demonstrates  the  feasibility  of  conducting 
simultaneous  NIRS  and  MRI  oximetry.  We  believe  the  correlation  of  tumour  vascular 
oxygenation  and  tumour  tissue  p02  can  provide  valuable  insights  into  tumour  pathophysiology 
and  response  to  interventions. 
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ABSTRACT 

1  The  goal  of  this  study  was  to  investigate  the  effect  of 
photothermal  laser  irradiation  on  rat  breast  tumor  (DMBA- 
4)  vascular  contents.  An  805-nm  diode  laser  was  used  in  our 
experiment  with  a  power  density  ranging  from  0.32  to  1.27 
W/cm2.  The  dynamic  changes  of  oxygenated  hemoglobin  and 
total  hemoglobin  concentrations,  A[Hb02]  and  A[Hb]total,  in 
rat  tumors  during  photothermal  irradiation  were  noninva- 
sively  monitored  by  a  near-infrared  spectroscopy  system.  A 
multichannel  thermal  detection  system  was  also  used  simul¬ 
taneously  to  record  temperatures  at  different  locations  within 
the  tumors.  Our  experimental  results  showed  that:  (1)  pho¬ 
toirradiation  did  have  the  ability  to  induce  hyperthermic  ef¬ 
fects  inside  the  rat  breast  tumors  in  a  single  exponential 
trend;  (2)  the  significant  changes  (P  <  0.005)  of  A[Hb02]  and 
A[Hb]total  in  response  to  a  low  dosage  of  laser  irradiation  (0.32 
W/cm2)  have  a  single  exponential  increasing  trend,  similar  to 
that  seen  in  the  tumor  interior  temperature;  and  (3)  the 
increase  in  magnitude  of  A[Hb02]  is  nearly  two  times  greater 
than  that  of  A[Hb]total,  suggesting  that  photoirradiation  may 
enhance  tumor  vascular  oxygenation.  The  last  observation 
may  be  important  to  reveal  the  hidden  mechanism  of 
photoirradiation  on  tumors,  leading  to  improvement  of  tumor 
treatment  efficiency. 


INTRODUCTION 

Laser  phototherapy,  a  nonsurgical  modality  for  cancer  treatment,  is 
gaining  widespread  acceptance  because  of  precise  energy  delivery 
into  the  tumor  tissue  (1-3).  The  laser  energy  delivered  to  the 
targeted  tumors  can  induce  localized  photomechanical,  photo¬ 
chemical  and  photothermal  reactions,  thus  killing  tumor  cells.  In 
general,  photochemical  reactions  may  cause  the  change  of 
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chemical  bonds  and  form  toxic  radicals,  such  as  release  of  singlet 
oxygen,  leading  to  the  death  of  organized  tissues.  Photomechanical 
reactions  may  induce  tissue  stress,  resulting  in  the  disruption  of 
tissue  cells  and  ejection  of  material.  Photothermal  reactions  may 
induce  hyperthermia  and  coagulation,  causing  cell  destruction  (4). 
In  particular,  thermal  cytotoxicity  effects  are  found  to  be  more 
profound  under  acidotic  (low  pH)  conditions,  which  are  often 
present  in  poorly  oxygenated  tissues  (5-7),  such  as  in  tumors. 
Many  research  studies  have  reported  the  existence  of  acute/chronic 
hypoxic  regions  within  the  majority  of  solid  tumors  (8-10).  As 
a  result,  cytotoxicity  induced  by  selective  photothermal  irradiation 
can  be  more  effective  in  tumor  tissue  than  in  normal  tissue.  To 
further  enhance  the  desirable  photothermal  effects,  selective 
absorbers  are  often  used  in  the  targeted  tumor  tissue  along  with 
laser  irradiation  to  cause  selective  and  localized  photon-tissue 
interactions  (11,12).  In  addition,  administration  of  an  immunoad- 
juvant  has  been  proved  to  increase  systemic  cancer  cure  and  long¬ 
term  resistance  to  cancer  of  the  same  origin,  when  combined  with 
the  selective  photochemical  treatment  (13,14).  Chen  et  al.  have 
developed  a  laser  immunotherapy  by  combining  laser  irradiation, 
a  laser-absorbing  dye  and  an  immunoadjuvant  to  improve  cancer 
therapy  efficiency  (1 1,12,15,16).  The  previous  experimental  results 
have  demonstrated  that  this  unique  approach  has  positive  effects  on 
both  primary  and  metastatic  tumors. 

Although  the  underlying  mechanism  of  the  laser  irradiation  has 
been  intensively  studied  (17-20),  the  hemodynamic  and  physio¬ 
logical  properties  inside  the  tumors  during  phototherapy  are  poorly 
understood.  Measuring  the  reduction  of  tumor  size  and  survival 
rates  is  the  current  practice  to  evaluate  the  therapeutic  outcome  of 
therapy.  Obviously,  noninvasive,  real-time  monitoring  of  tumor 
hemodynamic  characteristics  before,  during  and  after  therapy  is 
highly  desirable.  Such  a  technique  may  help  reveal  the  tumor 
physiology  and  therapeutic  process  caused  by  the  therapy, 
providing  treatment  prognoses  and  guidance  to  optimize  the  light 
dosimetry  so  as  to  improve  the  therapeutic  outcome.  Near-infrared 
spectroscopy  (NIRS)  has  been  demonstrated  in  our  recent  studies 
to  be  such  a  noninvasive  means  in  monitoring  tumor  vascular 
oxygenation  during  therapeutic  interventions  (21-23). 

In  this  study,  we  investigated  the  effects  of  photothermal  therapy 
on  tumor  vascular  contents  and  tumor  temperature.  An  805-nm 
diode  laser  was  used  for  photothermal  irradiation  of  mammary 
tumors  in  rats.  The  dynamic  changes  caused  by  the  laser  irradiation 
in  oxygenated,  deoxygenated  and  total  hemoglobin  concentrations, 
A[Hb02],  A[Hb],  and  A[Hb]total,  in  the  tumor  vasculature  were 
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Figure  1.  Experimental  setup  for  the  simultaneous  NIRS  and  thermal 
measurements  of  breast  tumors  under  photothermal  therapy.  PMT 
represents  a  photomultiplier  tube;  I&Q  is  an  in-phase  and  quadrature- 
phase  demodulator  for  retrieving  amplitude  and  phase  information.  The 
thermal  detection  system  comprises  multiple  thermal  sensing  probes,  three 
of  which  are  inserted  into  different  regions  of  the  tumor,  and  one  probe  was 
used  to  measure  room  temperature.  The  power-adjustable  laser  light  (805 
nm)  is  vertically  delivered  to  tumor  surface  with  a  2-cm  beam  diameter  for 
photoirradiation. 

monitored  by  an  NIRS  system.  A  multichannel,  thermal  monitor 
was  also  used  simultaneously  to  record  the  temperatures  at  three  dif¬ 
ferent  locations  inside  the  breast  tumors.  The  correlations  between 
the  optical  irradiation  dosages,  tumor  temperatures,  and  tumor  vas¬ 
cular  contents  were  studied.  The  dynamic  features  of  tumor  tem¬ 
perature  and  vascular  oxygenation  were  quantified  to  reveal  the 
dynamic  effects  within  the  tumors  under  the  laser  irradiation. 

MATERIALS  AND  METHODS 

NIRS.  The  changes  of  oxygenated,  deoxygenated  and  total  hemoglobin 
concentrations,  A[Hb02],  A[Hb]  and  A[Hb]totai,  in  tumor  tissue  caused  by 
laser  irradiation  were  determined  by  a  near-infrared,  homodyne,  frequency- 
domain  spectroscopy  system  (NIM,  Philadelphia,  PA),  which  has  been 
described  in  detail  elsewhere  (21-23).  Briefly,  the  amplitude-modulated 
light  at  140  MHz  at  wavelengths  of  758  nm  and  785  nm  illuminated  the 
tumor  surface  through  a  light-delivery  fiber  bundle  with  output  powers  of 
9  mW  and  11  mW  at  785  nm  and  758  nm  wavelengths,  respectively, 
measured  at  the  tip  of  the  bundle  by  an  optical  power  meter.  The 
experimental  setup  is  shown  in  Fig.  1.  The  805 -nm  laser  light  was  vertically 
applied  on  the  tumor  surface.  The  delivering  and  detecting  fiber  bundles  of 
the  NIRS  system  were  placed  horizontally  and  in  good  contact  with  the 
surface  of  the  tumors  in  a  transmittance  mode,  without  compressing  the 
tumors.  The  diffused  light  through  the  tumor  tissue  was  collected  and 
filtered  by  a  sharp  low-pass  optical  filter,  which  cut  off  the  optical  signal 
above  805  nm  from  the  laser  irradiation,  and  then  propagated  into 
a  photomultiplier  tube  (PMT)  through  the  detecting  fiber  bundle.  The  signal 
detected  by  the  PMT  was  demodulated  through  an  in-phase  and  quadrature- 
phase  (I&Q)  circuit,  and  the  amplitude  and  phase  of  the  signal  were 
recorded. 

Algorithms  for  calculations  of  A[HhO 2/,  A [Hh],  and  A[Hh]totai.  On  the 
basis  of  the  Beer-Lambert  law,  light  attenuation,  represented  by  optical 
density  (OD),  through  a  sample  is  proportional  to  the  concentration  of 
hemoglobin  inside  the  sample  (24-28), 

OD  =  Log(A0/A)  =  8  cl  (1) 

where  Aq  and  A  are  light  intensities  of  the  incident  and  transmitted  light, 
respectively,  £  is  the  extinction  coefficient  of  hemoglobin,  c  is  the 
concentration  of  hemoglobin,  and  /  is  the  length  of  light  path  through  the 
measured  sample.  When  the  measured  sample  has  a  mixture  of  oxygenated 
and  deoxygenated  hemoglobin,  Eq.  (1)  can  be  further  evolved  to  (28-30), 

OD(A.)  =  {£HbM[Hb]  +  8Hb02(^)[Hb02]}l  (2) 

where  OD(X)  is  the  optical  density  at  wavelength  X,  £Hb(X)  and  8Hb02(^)  are 
the  extinction  coefficients  at  wavelength  X  for  molar  concentrations  of 
deoxygenated  hemoglobin,  [Hb],  and  oxygenated  hemoglobin,  [Hb02], 
respectively,  assuming  ferrihemoglobin  is  minimal.  By  using  two  wave¬ 
lengths,  both  [Hb02]  and  [Hb]  can  be  determined  by  measuring  the  OD 


values  at  the  two  specific  wavelengths,  provided  that  the  values  for  £Hb(^) 
and  £Hb02(^)  are  known: 


^Hb (A,2)OD(A,i)  —  £Hb(Al  )OD(A,2) 

^[£Hb(^2)£Hb02(^l)  —  £Hb(^l) ^02(^2)]  ’ 

^Hb02  ( 'Xl )  OD  (  A,  j )  —  £Hb02(^“l)OD(X,2) 

^[£Hb(^l)£Hb02(^2)  —  £Hb(7<2)£Hb02(^l )] 


It  follows  that  changes  in  [Hb]  and  [Hb02]  can  be  consequently  given  as: 


pjbO  ]  —  £Hb(^2)AOD(^i)  —  £Hb(^l)AOD(?l2) 

^[£Hb(^2)£Hb02(^l)  —  ^Hb (^1  )£Hb02(^2)]  ’ 

_  £m>o2(A2)AOD(Xi)  —  sHb02(Ai)AOD(A,2) 

^[£Hb(^l)£Hb02(^2)  — 


(5) 

(6) 


where  AOD(X)  represents  a  change  in  optical  density  at  the  specific 
wavelength,  X,  and  equals  log^e/Ax).  AB  and  At  correspond  to  light 
intensities  measured  under  the  baseline  and  transient  conditions. 

Note  that  in  principle,  /  represents  the  optical  path  length  between  the 
source  and  detector.  Whereas  /  is  simply  the  physical  separation,  d,  between 
the  source  and  detector  through  a  nonscattering  medium,  exact  quantifi¬ 
cation  of  /  for  an  intact  tissue  or  organ  is  complex  because  of  light  scattering 
in  tissue.  Since  /  is  in  proportion  to  the  separation,  d,  we  can  associate  /  to 
d  as  /  =  DPF*d,  where  DPF  is  a  differential  path  length  factor  to  account 
for  light  scattering.  It  has  been  well  accepted  that  together  with  DPF, 
Eq.  (2)  can  be  treated  as  modified  Beer-Lambert  law;  and  consequently, 
Eqs.  (5)  and  (6)  can  be  correctly  used  to  quantify  changes  in  [Hb]  and 
[Hb02]  in  highly  scattering  media  (29,30),  such  as  in  intact  tissue  or  organs. 

To  be  consistent  with  our  previous  work,  we  adopted  in  this  paper  the  s 
values  published  by  Zijlstra  et  al.  (31).  We  had  to  interpolate  the  £  values  at 
the  two  wavelengths  used  in  our  study,  followed  by  certain  corrections  due 
to  the  interpolation  errors  by  phantom  calibration  measurements  (22). 
Specifically,  we  included  two  factors,  pi  and  p2,  for  the  calibrated 
algorithm,  as  given  below: 


A[Hb02] 


A[Hb] 


£Ht(7pi5mn)  X  OD(758  nm)  -  £Hb(758  nm)  X  OD(785  nm) 
/[Snb(785  nm)£Hbo2(758  nm)  -  £Hb(758  nm)£Hbo2(785  nm)]  ’ 

(7) 

sag nm>  X  OD(758  nm)  -  £Hbo2(758  nm)  X  OD(785  nm) 
/[£Hb(785  nm)£Hbo2(758  nm)  -  £Hb(758  nm)£Hbo2(785  nm)]  ’ 

(8) 


where  £Hb(758  nm)  =  1.418,  £Hb02(758  nm)  =  0.6372,  £Hb(785  nm)  = 
1.111,  and  £Hb02(785  nm)  =  0.766,  all  in  mM-1cm-1.  Note  that  a  factor  of 
four  has  been  multiplied  for  each  of  the  £s  at  the  respective  wavelengths  to 
account  for  light  absorption  from  four  hemes  per  hemoglobin  molecule  (34)  ?2 
since  the  extinction  coefficients  published  in  the  field  of  biochemistry  were 
expressed  on  a  heme  basis  (24-28).  Furthermore,  we  have  used  pt  =  1.103 
and  p2  =  0.9035  according  to  our  phantom  study  (22).  After  substituting  all 
of  the  parameters  into  Eqs.  (7)  and  (8),  we  have  arrived  at  the  final 
equations  to  quantify  changes  in  hemoglobin  concentration: 

a  riiur>  i  —2.658  •  OD(758  nm)  +  3.743  •  OD(785  nm) 

A[Hb02J  — - — - ,  (9) 

r  n  2.238  •  OD(758  nm)  -  1 .683  •  OD(785  nm)  ,  N 
A[Hb]  = - - - — - - - -.  (10) 


A[Hb]total  can  also  be  obtained  by  adding  Eqs.  (9)  and  (10), 

A[Hb]totll  =  A[Hb02]  +  A[Hb] 

_  — 0.42-OD(758  nm)  +  2.06-OD(785  nm) 


Notice  that  the  units  for  A[Hb02],  A[Hb],  and  A[Hb]total  are  mM; 
however,  the  values  obtained  from  Eqs.  (9)  to  (1 1)  are  scaled  by  a  factor  of 
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DPF.  Since  DPF  is  so  far  an  unknown  parameter  for  tumors,  we  include  it 
within  the  unit  as  mM/DPF. 

Laser  photothermal  irradiation  system.  The  system  used  for  the 
photothermal  therapy  consisted  of  an  NIR  diode  laser,  DIOMED  25 
(DIOMEDICS,  The  Woodlands,  TX),  with  an  emitting  peak  wavelength  at 
805  nm  and  a  maximum  power  output  of  25  W.  The  laser  light  was  coupled 
into  an  optical  fiber,  fitted  with  a  microlens  at  the  tip  (Pioneer  Optics, 
Windsor  Locks,  CT)  to  ensure  a  uniform  beam  density,  and  delivered  onto 
the  tumor  surface.  The  laser  beam  diameter  projected  on  the  tumor  surface 
was  set  as  2  cm  in  all  the  measurements  by  using  the  visible  aiming  light  in 
the  laser  system.  After  the  beam  diameter  was  adjusted,  the  aiming  light 
was  turned  off  during  NIR  measurements. 

Multichannel  thermal  monitor.  A  multichannel  thermal  monitor  (Omega, 
Stamford,  CT),  incorporating  an  OM-700  Omega  engineering  data 
acquisition  and  control  unit,  was  used  simultaneously  with  the  NIRS 
observation  to  monitor  the  internal  temperatures  at  several  locations  within 
the  tumors.  Specifically,  three  thermal  probes  were  inserted  at  the  top 
(immediately  below  the  overlaying  skin),  middle  and  bottom  positions  of 
the  tumors,  and  one  probe  was  located  outside  the  tumor  for  room 
temperature  recording.  The  data  of  each  probe  were  recorded  continuously, 
and  the  middle  thermal  probe  recorded  the  internal  tumor  temperature  near 
the  center  of  the  optical  field  of  NIRS. 

Animal  model.  The  transplantable,  metastatic  mammary  tumor  cells 
(DMBA-4)  (33-35)  were  implanted  in  one  of  the  inguinal  fat  pads  of 
female  Wistar  Furth  rats  (Harlan  Sprague  Dawley,  Indianapolis,  IN).  Once 
the  tumors  reached  0.5  to  0.8  cm  in  diameter,  the  rats  were  anesthetized  for 
the  experiments.  Hairs  around  tumors  were  removed  for  better  probe 
contact  in  the  NIRS  measurement  and  for  more  efficient  tumor  exposure  to 
the  laser  irradiation.  The  tumor  diameters  along  the  three  major  orthogonal 
axes  (a,  h,  c )  were  measured  to  determine  the  volume  of  the  tumors  by 
using  an  ellipsoid  approximation  with  the  formula  of  V  =  (n/6)abc. 

In  this  study,  16  rat  breast  tumors  were  used  in  two  different  treatment 
protocols.  Ten  rat  tumors  were  first  exposed  to  a  single  laser  irradiation  with 
a  laser  power  of  0.32  W/cm2  (1  W  laser  power  in  a  2-cm-diameter  beam) 
for  10  min.  Three  or  four  days  later,  four  of  the  treated  tumors,  along 
with  other  six  untreated  tumors,  were  exposed  to  repeated  laser  irradia¬ 
tions  with  a  power  density  of  0.32,  0.64,  0.96  and  1.27  W/cm2  (1  W,  2  W, 
3  W,  and  4  W  in  2-cm-diameter  beam)  for  10  min  with  a  recovery  period 
of  10  min  between  two  adjacent  irradiations.  Tumor  temperatures,  tumor 
vascular  A[HbC>2],  and  A[Hb]totai  were  measured  simultaneously  in  all 
the  experiments. 

Statistic  analysis.  A  Student’s  t-test  was  applied  in  data  analysis  to 
determine  whether  the  paired  samples  are  significantly  different.  We 
considered  P  <  0.05  as  the  criteria  for  significant  difference  between  the 
two  sample  groups.  In  addition,  a  mean  and  standard  deviation  (mean  ±  sd) 
were  usually  used  to  express  the  measured  data  and  their  dispersions. 

RESULTS 

Dynamic  responses  of  A[Hb02],  A[Hb]total 
and  temperature  to  photothermal  treatment 

Using  the  experimental  setup  shown  in  Fig.  1  and  Eqs.  (9)  and 
(1 1),  tumor  vascular  A[Hb02],  A[Hb]total  and  internal  temperatures 
were  quantified  simultaneously.  Figure  2  shows  the  temporal 
profiles  of  A[Hb02],  A [Hb]total  and  the  internal  temperatures  of 
a  rat  tumor  under  photothermal  irradiation  with  a  laser  power 
density  of  0.32  W/cm2.  The  experimental  data  consisted  of  four 
temporal  periods,  as  labeled  in  the  Fig.  In  period  1,  without 
irradiation,  the  thermal  probes  recorded  the  baseline  temperatures 
at  three  different  locations  inside  the  tumor  for  about  6  min.  In 
period  2,  the  NIRS  system  was  turned  on  while  the  thermal 
recording  continued,  and  the  recorded  thermal  readings  showed  no 
obvious  perturbation  from  the  NIRS.  After  the  NIRS  signals  were 
stabilized,  a  10-min  laser  irradiation  started.  During  the  laser 
treatment,  the  tumor  temperatures,  A[Hb02],  and  A[Hb]total 
showed  significant  elevation  (period  3  in  Fig.  2a).  Tumor 
temperatures  increased  significantly  ( P  <  0.0001)  with  the 
irradiation.  The  tumor  volumes  near  the  tumor  surface  had  more 
significant  thermal  effects:  the  temperature  at  the  top  (T-top)  varied 


from  25.3  ±  0.1°C  to  30.3  ±  0.1°C  (mean  ±  sd),  while  the 
temperatures  at  the  middle  and  bottom  (T-middle  and  T-bottom)  of 
the  tumor  were  elevated  from  26.1  ±  0.1°C  to  29.8  ±  0.1°C  and 
from  27.5  ±  0.1°C  to  28.9  ±  0.1°C,  respectively.  The  maximal 
changes  in  T-top,  T-middle  and  T-bottom  caused  by  the  laser 
irradiation  were  5.0  ±  0.1°C,  3.7  ±  0.1°C  and  1.4  ±  0.1°C, 
respectively.  Similarly,  both  A[Hb02]  and  A[Hb]total  displayed 
significant  increases  over  the  entire  irradiation  period  (P  <  0.005), 
with  the  mean  elevations  of  0.04  ±  0.01  mM/DPF  and  0.02  ± 
0.01  mM/DPF,  respectively.  Notice  that  the  maximal  values  of 
A[Hb02]  and  A[Hb]total  shown  in  Fig.  2a  are  ~0.1  and  0.05  mM/ 
DPF.  When  the  laser  beam  was  turned  off  in  period  4,  A[Hb02], 
A[Hb]total  and  temperatures  all  decreased,  with  the  temperatures 
reaching  the  baselines  quickly  and  A[Hb02]  and  A[Hb]total 
declining  more  gradually  and  not  quite  reaching  the  baselines 
within  the  10-min  recovery  period. 

To  quantify  the  dynamic  behaviors  observed  in  Fig.  2a,  single¬ 
exponential  functions  were  applied  to  T-top,  T-middle,  A[Hb02], 
and  A[Hb]totai  on  both  rising  and  falling  periods,  as  shown  in  Fig. 
2b, c.  The  time  constants  obtained  from  such  fittings  reveal  the 
dynamic  responses  of  signals  to  the  initiation  and  termination  of 
photothermal  irradiation.  As  shown  in  Fig.  2b,  T-top  (time  constant 
of  4.5  ±  0.5  min)  has  a  faster  response  to  laser  irradiation  than 
T-middle  (time  constant  as  14.4  ±  0.6  min).  The  time  constants  for 
A[Hb02]  and  A[Hb]total  are  3.85  ±  0.03  min  and  3.95  ±  0.06  min, 
respectively.  In  response  to  the  termination  of  laser  irradiation, 
T-top  and  T-middle  have  a  time  constant  of  3.9  ±  0.5  min  and 
12.3  ±  0.6  min,  respectively  (Fig.  2c).  It  is  clearly  seen  that  the 
dynamic  thermal  responses  of  the  tumor  to  laser  irradiation  and  to 
the  termination  of  laser  irradiation  are  well  matched  both  at  the 
tumor  surface  and  within  the  tumor.  However,  the  time  constants 
for  A[Hb02]  and  A[Hb]total  to  return  to  the  baseline  are  much 
longer  in  the  falling  period,  with  18.5  ±  0.5  min  and  24.3  ± 
0.8  min,  respectively,  also  shown  in  Fig.  2c. 

Similar  protocols  and  measurements  were  performed  on  the 
other  nine  tumors.  Figure  3  shows  the  average  increases  of  the 
temperature  at  different  locations.  Temperatures  near  the  tumor 
surface  have  maximal  increases,  whereas  local  temperatures  near 
the  bottom  of  the  tumors  have  minimal  increases  observed  from  all 
the  10  tumors.  The  average  elevations  for  T-top  and  T-bottom  are 
7  ±  3°C  and  3  ±  2°C,  respectively.  The  thermal  readings  from  the 
middle  position  (T-middle)  represent  the  average  temperatures  in 
the  tumor  volumes  probed  by  the  NIRS  beams,  and  the  average 
increase  in  T-middle  for  the  10  rats  is  5  ±  2°C.  The  statistical 
analysis  indicates  significant  differences  between  T-top,  T-middle, 
and  T-bottom  (P  <  0.05). 

To  investigate  dynamic  behaviors  of  the  tumors  in  response  to 
the  initiation  and  termination  of  laser  irradiation,  the  time  constants 
of  T-top,  T-middle,  A[Hb02],  and  A[Hb]total  obtained  from  the 
rising  and  falling  periods  are  summarized  in  Table  1.  The  results 
demonstrate  several  points:  (1)  an  average  dynamic  change  in 
T-top  is  significantly  (P  <  0.05)  faster  than  that  in  T-middle  in  the 
rising  period,  with  a  time  constant  of  3.7  ±  1.1  min  for  T-top  and 
6.4  ±3.2  min  for  T-middle.  A  similar  significant  difference  (P  < 
0.01)  also  holds  in  the  falling  period  with  a  time  constant  of  3.8  ± 
1.2  min  for  T-top  and  6.6  ±  2.7  min  for  T-middle,  respectively.  (2) 
The  dynamic  responses  of  A[Hb02],  A[Hb]total,  and  T-middle  have 
no  significant  differences  during  the  initiation  of  photoirradiation 
(P  >  0.4).  (3)  The  average  time  constant  of  T-top  in  the  rising 
period  is  similar  (P  >  0.9)  to  that  in  the  falling  period,  and  so  is 
that  of  T-middle  (P  >  0.7).  However,  (4)  the  dynamic  response  of 
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Figure  2.  (a)  Temporal  profiles  of  the  changes  in  tumor  oxygenation  and 
total  hemoglobin  concentrations,  A[Hb02]  and  A[Hb]totai,  and  three 
intratumor  temperatures  during  photothermal  intervention  in  a  representative 
DMBA-4  rat  breast  tumor  (2.4  cm3).  The  profile  is  divided  into  four 
periods,  and  the  laser  irradiation  was  applied  only  in  period  3.  (b) 
Respective  single-exponential  curve  fittings  for  the  rising  period  of  T-top, 
T-middle,  A[Hb02]  and  A[Hb]tota|  displayed  in  (a).  The  top  .v-axis  is 
used  for  the  T-top  and  T-middle  curves,  whereas  the  bottom  .v-axis  is 
for  A[Hb02]  and  A[Hb],ota|.  Single  exponential  fitting  yields:  T-top  = 
5.76*(1  -  exp[— {f  -  147}/4.49])  +  24.95  with  R  =  0.99;  T-middle  = 


T-top 


T-middle 


T-bottom 


Positions  in  tumors 

Figure  3.  Average  changes  in  temperature  from  the  10  rats  for  T-top, 
T-middle  and  T-bottom. 


A[Hb02]  in  the  falling  period  (a  time  constant  of  28  ±  22  min)  is 
significantly  slower  ( P  <  0.01)  than  that  in  the  rising  period  (a  time 
constant  of  5.3  ±  2.8  min),  and  so  is  A[Hb]total  response,  with 
a  time  constant  of  6.2  ±  3.7  min  and  41  ±  32  min  for  the  rising 
and  falling  periods,  respectively  (P  <  0.01). 

Furthermore,  although  A[Hb02]  and  A[Flb]total  have  similar 
dynamic  trends  in  response  to  photoirradiation,  a  difference  be¬ 
tween  their  magnitudes  exists.  A  strong  correlation  (R2  =  0.95)  be¬ 
tween  A[Flb02]  and  A[Hb]totai  has  been  observed,  as  shown  in  Fig. 
4,  where  the  amplitude  in  A[Hb02]  is  nearly  twice  that  of  A[Hb]total 
for  all  10  tumors  (with  a  slope  of  the  fitted  curve  being  0.54). 


Effect  of  the  optical  power  density  on  tumor 
A[Ffb02],  A[Hb]total  and  temperatures 

Repeated  laser  irradiation  was  applied  to  the  10  tumors,  4  of  which 
were  used  in  the  previous  single-irradiation  treatment  3  or  4  days 
earlier.  Figure  5a  shows  the  profiles  of  the  tumor  vascular 
A[Hb02],  A[Hb]totai  and  tumor  temperatures  in  response  to 
irradiations  under  power  densities  of  0.32,  0.64,  0.96  and  1.27 
W/cm2,  with  an  irradiation  period  of  10  min  and  a  recovery  period 
of  10  min  between  each  treatment.  Tumor  temperatures  (T-top, 
T-middle  and  T-bottom)  increased  significantly  (P  <  0.0001)  from 
the  baseline  to  different  maximum  values  under  different  laser 
powers.  The  maximal  changes  in  temperature  in  the  middle  of  the 
tumor  (T-middle)  due  to  consecutive  irradiation  displayed  a  strong 
linear  dependence  on  the  laser  power  density,  as  shown  in  Fig.  5b. 
However,  A[Hb02]  and  A[Hb]total  showed  interesting  behaviors  in 
response  to  the  thermal  irradiation  at  different  light  dosages.  Under 
relatively  low  laser  power  irradiations  (0.32,  0.64,  0.96  W/cm2), 
both  A[Hb02]  and  A[Hb]totai  increased  significantly  (P  <  0.0001) 
with  respect  to  the  baselines,  and  they  both  declined  significantly 
under  a  high  laser  power  irradiation  (1.27  W/cm2),  as  shown  in 


7.28*(1  -  exp[-{f  -  14}/14.3])  +  25.98  with  R  =  0.99:  A[HbOJ  = 
0.093*(1  -  exp[-{t  -  14}/3.85])  +  0.013  with  R  =  0.98;  A[Hb]total  = 
0.044*(1  —  exp[— {t  —  14}/3.95])  +  0.012  with  R  =  0.98;  (c)  Respective 
single-exponential  curve  fittings  for  the  falling  period  of  T-top,  T-middle, 
A[Hb02]  and  A[Hb]totai  displayed  in  (a).  Single  exponential  fitting 
yields:  T-top  =  -4.48*(1  -  exp[-{f  -  25}/3.86])  +  29.23  with  R  =  0.98"; 
T-middle  =  -3.83*(1  -  exp[-{f  -  25}/12.29])  +  30.42  with  R  =  0.99; 
A[Hb02]  =  -0.163*(1  -  exp[-{t  -  25J/18.5])  +  0.10  with  R  =  0.94; 
A[Hb]total]  =  -0.1 14*(1  -  exp[— {f  -  25}/24.3])  +  0.06  with  R  =  0.92. 


Photochemistry  and  Photobiology,  2005,  81 


Table  1.  Time  constants  (minutes)  of  T-top,  Tmiddle.  A[Hb02],  and  A[Hb]total  with  respect  to  photoirradiation  and  tennination 


Tumor 

size(cm3) 

Rising  part  (in  response  to  irradiation) 

Falling  part  (in  response  to  termination) 

T-top 

Tmiddle 

A[Hb02] 

A[Hb]total 

T-top 

Tmiddle 

A[HbO,] 

A[Hb]total 

2.4 

4.49 

14.3 

3.85 

3.95 

3.86 

12.29 

18.5 

24.3 

2.5 

2.47 

6.77 

1.99 

1.53 

3.17 

8.27 

50 

74 

0.68 

4.83 

4.83 

1.76 

1.66 

3.87 

4.45 

20 

25 

0.53 

4.61 

7.62 

1.84 

1.91 

4.56 

8.37 

3.63 

2.75 

1.2 

3.28 

5.17 

6.9 

9.8 

4.20 

6.30 

80 

120 

3.3 

4.70 

6.65 

7.34 

7.83 

3.58 

4.45 

11.3 

29.9 

1.9 

1.89 

3.85 

10.0 

12.0 

1.30 

3.43 

28 

37 

3.4 

2.69 

3.13 

6.76 

7.81 

2.96 

4.58 

35 

45 

2.1 

3.25 

4.31 

7.32 

7.58 

4.52 

5.65 

23.5 

30 

3.0 

4.50 

7.80 

5.60 

8.40 

5.50 

8.30 

15.5 

23.4 

Average 

3.7 

6.4 

5.3 

6.2 

3.8 

6.6 

28. 

41 

SD* 

1.1 

3.2 

2.8 

3.7 

1.2 

2.7 

22. 

33 

*SD  =  standard  deviation. 


Fig.  5b.  Similar  to  the  feature  observed  in  single  irradiation, 
A[Hb02]  values  had  greater  magnitude  than  A[Hb]total,  whereas 
their  dynamic  trends  were  consistent. 

To  investigate  the  accumulative  effects  on  A[Hb02]  and 
A[Hb]total  from  several  thermal  irradiations  given  10  min  apart 
(Fig.  5a),  we  had  also  taken  a  direct  reading  for  A[Hb02]  and 
A[Hb]total  with  a  high  laser  power  of  1.27  W/cm2  without  any 
preirradiation.  Figure  5c  shows  an  example  of  this  case,  demon¬ 
strating  that  A[Flb02]  indeed  decreases  after  an  initial  increase  for 
a  few  minutes,  whereas  A[Flb]total  has  less  trend  to  drop  within 
the  10-min  treatment.  The  similarity  for  the  ratio  of  A[Hb02]  to 
A[Hb]totaj  close  to  a  factor  of  2  still  exists. 

Temperature  data  for  all  10  tumors  under  different  power 
irradiations  are  plotted  in  Fig.  6a.  The  maximum  changes  in  T- 
rniddle  increase  linearly  with  the  optical  radiation  dosages,  with 
a  correlation  of  R2  >  0.91  for  all  the  cases  and  having  an  average 
slope  of  13  ±  3  (°C/[W/cm2]).  To  distinguish  the  tumors  pretreated 
and  nonpretreated  with  a  single  low-power  irradiation  (0.32 
W/cm2)  3  or  4  days  before,  the  slopes  of  the  linear  regression  for 
each  of  the  tumors  shown  in  Fig.  6a  were  classified  into  two 
groups,  i.e.  without  pretreatment  and  with  pretreatment  (Fig.  6b). 
A  Student’s  ?-test  showed  that  there  is  no  significant  difference 
( P  =  0.95)  between  the  two  groups  of  slopes. 

DISCUSSION 

In  this  study,  we  have  investigated  vascular  A[Hb02],  A[Hb]total 
and  internal  temperature  in  response  to  photothermal  therapy 
simultaneously  in  rat  breast  tumors  by  using  an  NIRS  system  and 
a  thermal  monitor,  respectively.  The  results  obtained  in  this  study 
clearly  demonstrate  the  compatibility  of  the  NIRS  system  with  the 
temperature  monitoring  system,  without  noticeable  interference 
between  the  two  systems  (Fig.  2).  The  NIRS  system  offers  a  real¬ 
time  measurement  of  A[Hb02]  and  A[Hb]totai  in  tumor  vasculature 
within  the  optical  field  of  the  NIR  probes,  while  the  multichannel 
thermal  monitor  records  the  temperatures  at  specific  locations 
within  the  tumors  in  real  time.  The  compatibility  of  the  two 
systems  permits  simultaneous  determinations  of  thermal  and 
vascular  characteristics  of  treated  tumors  under  photothermal 
therapy,  providing  valuable  insight  into  the  dynamic  relationship 
between  the  tumor  vascular  contents  and  the  thermal  effect  in 
treated  tumor  tissues. 


Thermal  treatment  to  tumors  has  been  demonstrated  to  cause 
cellular  cytotoxicity,  and  thus,  to  destroy  the  tumor  cells  (36-38). 
Laser  irradiation  applied  in  this  study  proves  its  ability  to  induce 
photothermal  effects  inside  the  rat  breast  tumors.  The  single 
exponential  increase  in  tumor  temperature  obtained  in  our  results 
indicates  that  an  extended  exposure  period  to  the  laser  beam  may 
gradually  stabilize  the  tumor  temperature  at  a  certain  level, 
depending  on  both  the  laser  power  density  and  the  individual 
tumors  (as  seen  in  Figs.  2a  and  5a).  Manipulation  of  the  laser 
power  may  help  adjust  the  internal  tumor  temperature,  resulting  in 
optimal  tissue  destruction.  Although  a  strong  linear  relation 
between  the  maximum  temperature  increase  within  the  tumors 
and  the  irradiation  power  density  has  been  observed  for  all  the 
tumors,  the  actual  dependence  of  tumor  temperature  on  the  laser 
power  density  does  vary  from  tumor  to  tumor  (Fig.  6).  Such 
intertumor  variability  demonstrates  the  necessity  of  real-time 
monitoring  for  tumor  temperature  under  photothermal  therapy. 
Moreover,  the  tumor  temperatures  near  the  surface  (T-top)  respond 
to  the  photoirradiation  more  rapidly  and  with  a  greater  increase 
than  those  within  deeper  tumor  tissue  (T-middle)  (as  seen  in  Figs. 
2,  3.  and  5a,  and  Table  1).  Such  large  thermal  effects  can  cause 
tissue  damage  at  the  tumor  surface.  To  avoid  such  thermal  damage. 


Figure  4.  Relationship  between  AfHbOJ  and  A[Hb]totai  for  the  10  rat 
tumors  under  single  laser  irradiation  treatment  (0.32  W/cm2  for  10  min) 
with  a  linear  fit  of  A[Hb]total  =  0.54*A[HbO2]  —  0.0039  (the  solid  curve) 
and  a  correlation  of  R2  =  0.96.  The  dashed  line  represents  the  situation  of 
A[Hb02]  =  A[Hb],otal. 


Yueqing  Gu  et  al. 


without  pretreatment  with  pretreatment 

Figure  6.  (a)  Correlations  between  the  maximum  changes  in  T-middle 
reached  at  the  end  of  each  photothermal  treatment  cycle  and  the  optical 
irradiation  power  density  from  the  10  treated  rats.  The  numbers  after  the 
symbol  “#”  are  labels  for  different  tumors.  The  numbers  inside  the 
parentheses  indicate  the  tumor  volume,  with  unit  of  cm3.  Note  that  there  are 
four  tumors  that  had  preirradiations  3—4  days  earlier,  and  their  tumor 
volumes  have  been  possibly  changed  in  comparison  to  those  listed  in 
Table  1.  (b).  The  slopes  in  (a)  for  all  the  tumors  were  classified  into  two 
groups:  without  pretreatment  and  with  pretreatment.  A  f-test  between  these 
two  groups  of  slopes  gives  rise  to  P  =  0.95. 


Time(mln) 

Figure  5.  (a)  The  temporal  profiles  of  vascular  A[Hb02].  A[Hb]tota|  and 
temperatures  of  a  tumor  (#6,  3.0  cm3)  in  response  to  repeated  laser 
irradiations  with  different  laser  power  densities  (0.32,  0.64,  0.96  and  1.27 
W/cm2)  for  a  10-min  treatment  cycle  with  a  10-min  recovery  period 
between  each  treatment,  (b).  The  correlations  between  the  changes  in 
A[Hb02]  (solid  triangles).  A[Hb]totai  (squares),  T-middle  (open  triangles), 
and  laser  irradiation  power  densities  in  the  tumors.  The  thick  straight  line  is 
the  best-fitting  curve  for  the  tumor  temperature  increases  as  function  of 
laser  dosage,  (c)  The  temporal  profiles  of  vascular  A[Hb02]  and  A[Hb]totai 
of  a  tumor  irradiated  with  a  single  laser  power  density  of  1.27  W/cm2. 

intratumoral  injection  of  selective  absorbers,  such  as  indocyanine 
green,  has  been  used  by  Chen  et  al.  to  selectively  manipulate  the 
temperature  within  the  tumor  without  severe  damage  at  the  tumor 
surface  (1 1,12,39). 

It  is  seen  clearly  in  this  study  that  tumor  vascular  contents,  i.e. 
oxygenated  and  deoxygenated  hemoglobin,  are  greatly  affected  by 
photothermal  therapy.  When  the  photothermal  irradiation  with 
a  low  power  density  (0.32  W/cm2)  was  applied  to  the  breast 
tumors,  both  A[HbC>2]  and  A[Hb]totai  increased  significantly  (P  < 
0.005)  in  all  the  tumors  (Figs.  2a,  4  and  5a).  The  dynamic 
behaviors  of  A[Hb02]  and  A[Hb]totai  in  response  to  photo¬ 
irradiation,  quantified  by  the  time  constants,  are  in  the  same  order 
as  those  of  tumor  interior  temperature,  T-middle,  while  their 


behaviors  in  response  to  the  termination  of  photoirradiation  are 
much  slower  than  those  of  T-middle  (Fig.  2  and  Table  1).  In 
contrast  to  the  lower  laser  power  irradiation,  higher  laser  power 
densities  could  cause  interesting  responses  in  tumor  vascular 
contents.  As  seen  in  Fig.  5b, c,  A[Hb02]  started  to  decrease  when 
the  laser  power  approached  a  certain  level,  either  with  or  without 
preirradiation.  Intriguingly,  although  A[Hb02]  has  the  same 
dynamic  trend  as  A[Hb]total  in  response  to  photoirradiation,  the 
magnitude  of  A[Hb02]  is  more  profound  than  that  of  A[Hb]total 
(Figs.  4,  5). 

The  analysis  on  the  dynamic  behaviors  of  tumor  vascular  contents 
and  interior  tumor  temperature  in  response  to  photoirradiation 
allows  us  to  better  understand  the  mechanism  of  photothermal 
therapy.  The  time  constants  of  A[Hb02]  as  well  as  of  A[Hb]total 
have  the  same  order  as  those  of  interior  tumor  temperature.  We  may 
speculate  that  an  increase  in  tumor  temperature  may  lead  to  vessel 
dilation  within  tumor  vasculature,  and  thus,  resulting  in  an  increase 
in  blood  inflow.  Such  an  increase  in  blood  inflow  will  give  rise 
to  an  increase  in  A[Hb]total  as  well  as  A[Hb02]. 

Besides  thermal  effects,  it  is  reasonable  to  expect  other 
photoactivation  processes  involved  in  the  laser  irradiation.  It  is 
known  that  direct  exposure  of  vascular  tissue  to  visible  light 
causes  vasodilation  (40),  whereas  the  thermal  effect  from  visible 
light  is  minor.  Photoactivation  has  also  been  reported  to  cause 
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upregulation  of  nitric  oxide,  and  thus,  produce  vasodilation  (41 
43).  The  increase  of  tumor  blood  contents  due  to  laser  irradiation 
observed  in  our  study  may  be  attributed  to  the  vessel  dilation 
induced  by  the  direct  effect  of  heating  as  well  as  photoactivation. 
Consistently,  the  single-exponential  trend  of  both  A[Hb]total  and 
A[Hb02]  in  response  to  laser  irradiation  indicated  the  existence 
of  vessel  dilation.  The  fact  that  the  time  constants  of  A[Hb02]  and 
A[Hb]total  in  response  to  the  laser  termination  are  longer  than  those 
in  response  to  the  laser  initiation  (Table  1,  Figs.  2a,  5a)  implies 
a  longer  settling  time  needed  for  the  recovery  of  dilated  vessels. 

In  the  laser  thermal  treatment,  we  realized  that  A[Hb02]  as  well 
as  A[Hb]totai  did  not  completely  return  to  the  baseline  during  each 
of  the  recovery  periods  (Figs.  2a,  5a).  Figure  5c  was  taken  from 
a  tumor  that  was  exposed  to  a  single  high-power  laser  irradiation 
(1.27  W/cm2),  displaying  somewhat  different  behavior  in  compar¬ 
ison  to  that  under  consecutive  repeated  treatments:  the  A[Hb02] 
signals  increased,  followed  by  an  obvious  decrease.  These  phe¬ 
nomena  lead  us  to  a  speculation  that  the  previous  laser  exposure 
“preconditions”  the  treated  tumor,  with  respect  to  consecutive 
treatments,  and  affects  the  response  of  tumor  to  subsequent  high 
irradiation  exposure,  although  the  single  lower-power  irradiation 
pretreated  3  or  4  days  earlier  has  no  obvious  influence  on  the  later 
treatment  (Fig.  6b).  It  is  reasonable  to  believe  that  the  “precondi¬ 
tioning”  can  be  minimal  if  the  time  intervals  between  two  repeated 
irradiations  are  two  to  three  times  longer  than  the  time  constants 
during  the  recovery  periods,  such  as  in  the  order  of  1—1.5  h. 

The  interesting  phenomenon,  i.e.  changes  of  oxygenated 
hemoglobin  concentration  are  nearly  two  times  greater  than  those 
of  total  hemoglobin  concentration,  allures  us  to  speculate  that  the 
increase  in  A[Hb02]  comes  partly  from  the  increase  of  A[Hb]totai, 
and  partly  from  the  conversion  of  deoxygenated  hemoglobin  to 
oxygenated  hemoglobin.  The  correlation  of  magnitude  between 
A[Hb02]  and  A[Hb]total  unambiguously  demonstrates  that  the 
increased  total  hemoglobin  concentration,  due  to  blood  vessel 
dilation,  is  mostly  oxygenated  after  the  irradiation  starts.  It 
basically  implies  that  photoirradiation  may  increase  the  tumor 
oxygenation,  which  will  be  beneficial  to  some  nonsurgical  treat¬ 
ments,  such  as  chemotherapy  and  radiotherapy  (37). 

In  summary,  changes  in  breast  tumor  temperature  and  vascular 
oxygenation  have  been  simultaneously  measured  using  a  multi¬ 
channel  thermal  monitor  and  an  NIRS  system,  while  the  tumors 
were  under  photothermal  irradiation.  The  results  have  demonstrat¬ 
ed  that:  (1)  the  compatibility  of  thermal  monitoring  system  and  the 
NIR  system  permits  simultaneous  determinations  of  thermal  and 
vascular  characteristics  of  the  treated  tumors;  (2)  photoirradiation 
did  have  the  ability  to  induce  thermal  effects  inside  the  rat  breast 
tumors  in  a  single  exponential  trend,  thus  leading  to  destruction  of 
the  tumor  cells;  (3)  the  response  of  A[Hb02]  and  A[Hb]totai  to 
photoirradiation  may  be  attributed  to  vessel  dilation  under  both 
thermal  and  photo  effect;  (4)  low-power  laser  irradiations 
predelivered  3  or  4  days  earlier  have  no  effect  on  the  later 
treatments;  the  previous  laser  exposures  within  a  time  interval  of 
10  min  do  influence  the  response  of  tumor  vascular  contents  to 
subsequent  high  irradiation;  (5)  the  fact  that  the  change  in  the 
magnitude  of  A[Hb02]  is  nearly  two  times  greater  than  that  of 
A[Hb]total  suggests  that  photoirradiation  may  enhance  tumor 
vascular  oxygenation.  The  primary  results  from  this  study  clearly 
indicate  that  besides  a  considerable  thermal  effect,  low-power 
density  irradiation  could  also  result  in  a  significant  enhancement  of 
tumor  vascular  oxygenation;  this  enhancement  may  lead  to 
improvement  of  treatment  efficiency. 
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Extinction  coefficients  of  hemo¬ 
globin  have  been  studied  for  five 
decades  by  clinical  chemists  and 
biochemists,  particularly  for  lab¬ 
oratory  spectrophotometric  measure¬ 
ments.  A  representative  list  of  literature 
can  be  found  in  [1] — [12],  and  the  studies 
on  the  temperature  dependence  of 
extinction  coefficients  of  hemoglobin 
can  be  found  in  [13]  and  [14]. 

However,  during  the  last  ten  to  15 
years,  near  infrared  spectroscopy  (NIRS) 
and  imaging  for  tissue  vascular  oxygena¬ 
tion,  breast  tumor  detection,  and  func¬ 
tional  brain  imaging  have  been 
intensively  developed  for  in  vivo 
measurements  by  groups  of  physicists, 
biomedical  engineers,  and  mathemati¬ 
cians.  In  the  approach  of  NIRS,  NIR 
light  in  the  wavelength  range  of  650-900 
nm  is  utilized  to  illuminate  tissue  in  vivo, 
and  the  transmitted  or  reflected  light 
through  tissue  is  recorded  for  the  quan¬ 
tification  of  hemoglobin  concentrations 
of  the  measured  tissue  vasculature.  In 
order  to  achieve  mathematical  conver¬ 
sion  from  the  detected  light  intensity  at 
different  wavelengths  to  hemoglobin 
concentration,  extinction  coefficients  of 
hemoglobin,  e,  must  be  used. 

While  the  engineers  and  physicists 
working  in  the  NIR  field  have  found  the 
correct  e  values  to  use,  there  has  been 
controversy  on  what  e  values  should  be 
used  for  in  vivo  NIRS  in  comparison 
with  the  conventional  e  that  most  bio¬ 
chemists  have  used  in  the  laboratories 
for  in  vitro  measurements.  The  purpose 
of  this  article  is  to  address  this  issue  and 
help  biomedical  engineers  and  physi¬ 
cists  gain  a  better  understanding  of  e  to 
be  used  for  NIRS  and  NIR  imaging. 

Spectrophotometric 
Measurements  to  Determine 
Hemoglobin  Concentration 

Hemoglobin  is  a  molecule  in  the  red 
blood  cells  that  has  a  role  of  deliver¬ 


ing  oxygen  to  tissue  cells. 
Hemoglobin  is  composed  of  four 
heme  groups  and  a  protein  group, 
known  as  a  globin.  Historically,  for 
spectrophotometric  experiments,  bio¬ 
logical  chemists  and  biochemists  uti¬ 
lized  Beer-Lambert’s  law  and 
developed  the  notation  of  absorbance 
to  express  light  absorption  as  a  func¬ 
tion  of  hemoglobin  concentration  as 
given  in  [  1  ]— [5] : 

OD  =  Log  (Iq/I)  =  sc  L  (1) 

where  OD  is  the  optical  density,  Ig  is 
the  light  intensity  of  incident  light,  /  is 
the  light  intensity  of  transmitted  light,  s 
is  the  extinction  coefficient  of  hemoblo- 
gin,  c  is  the  concentration  of  hemoglo¬ 
bin,  and  L  is  the  length  of  light  path 
through  solution. 

When  the  measured  sample  has  a 
mixture  of  oxygenated  and  deoxy- 
genated  hemoglobin,  (1)  can  be  further 
expanded  as  [3] — [5], 

ODx  =  [£Hb"[Hb]  +  £Hbo/[Hb02]}L 

(2) 

where  OD^  is  the  optical  density  or 
absorbance  at  wavelength  X  and 
fiHbW  and  £m>o2W  are  the  extinction 
coefficients  at  wavelength  X  for 
molar  concentrations  of  deoxygenat- 
ed  hemoglobin,  [Hb],  and  oxygenated 
hemoglobin,  [Hb02],  respectively, 
assuming  ferrihemoglobin  is  mini¬ 
mal.  In  some  references,  the  light 
path,  L,  was  taken  as  1  cm  without 
mentioning  the  definition  of  L  each 
time  [3] — [5] .  Both  [Hb02]  and  [Hb] 
can  be  determined  by  measuring  the 
light  absorbance  at  the  two  specific 
wavelengths,  provided  that  the  values 
for  £Hb(X)  and  £Hbo2M  are  known,  as 
expressed  below. 


[Hb02] 

£HbX2ODM  ~  £HbA'ODA2 

LLfiHb^eHbOj^1  —  fiffi/'fiHbO,12) 

(3) 

[Hb] 

£Hbo2A2ODXl  -  £Hbo2A'ODA2 

L(£ffl,A'1£Hb02;‘2  —  fiHb^fiHbO,11) 

(4) 

It  follows  that  changes  in  [Hb]  and 
[Hb02]  can  be  consequently  given  as 

A[Hb02] 

£HbA2AODM  -£HbMAODi2 
LtfiHb^eHbO;/1  —  £HbX2£Hb02l2) 

(5) 

A[Hb] 

_  £Hbo2X2AODA|  -  £Hbo2X|  AODX2 

L(£HbM£Hbo2A2  -  £Hb,'2£Hbo2Al) 

(6) 


A[Hb]total  =  A[Hb]  +  A[Hb02],  (7) 

where  AODA  represents  a  change  in 
optical  density  at  the  specific  wave¬ 
length,  X,  and  equals  logO^/Ix).  Ib  and 
IT  correspond  to  light  intensities  mea¬ 
sured  under  the  baseline  and  transient 
conditions. 

These  spectrophotometric  calcula¬ 
tions  seem  straightforward  mathe¬ 
matically  and  have  been  used  for 
several  decades  by  biochemists  to 
quantify  [Hb]  and  [Hb02]  in  the  lab¬ 
oratory  measurements.  However, 
close  attention  must  be  paid  to  the 
definition  and  accuracy  of  s  since 
biochemical  details  in  obtaining  £ 
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values  give  rise  to  the  different  quan¬ 
tification  of  hemoglobin  concentra¬ 
tion.  In  early  publications  on 
spectrophotometry  of  hemoglobin 
[3],  it  was  clearly  stated  that  all 
extinction  coefficients  were 
expressed  on  a  heme  basis,  using  a 
term  of  equivalent,  where  per  equiv¬ 
alent  of  hemoglobin  was  assumed  to 
be  66,800/4  or  16,700  gm,  i.e.,  one- 
quarter  of  the  molecular  weight  of 
the  hemoglobin  molecule  [12].  While 
spectroscopic  absorption  measure¬ 
ments  became  a  popular  methodolo¬ 
gy  for  biochemists  to  quantify  [Hb] 
and  [Hb02]  in  laboratories  [5],  [6], 
[12]  the  notation  of  extinction  coeffi¬ 
cients  being  based  on  one  heme 
group  (or  per  equivalent)  gradually 
has  not  been  mentioned  and  has 
become  conventionally  understood 
by  biological  chemists  for  the  last 
two  to  three  decades. 

Near-Infrared  Spectroscopy  Used 
to  Determine  Hemoglobin 
Concentration 

For  the  last  15  or  more  years,  signifi¬ 
cant  research  efforts  have  been  con¬ 
ducted  using  NIRS  for  in  vivo 
quantification  of  hemoglobin  concen¬ 
trations  from  intact  living  tissues.  At 
the  early  stage  of  development,  direct 
utilization  of  (3)  and  (4)  was  per¬ 
formed  to  obtain  values  or  changes  of 
oxygenated  and  deoxygenated  hemo¬ 
globin  concentration  from  living  tis¬ 
sues  in  vivo  [15]— [18],  without 
particularly  noticing  that  most  of  the  e 
values  published  in  biological  litera¬ 
ture  were  extinction  coefficients  for 
hemoglobin  per  equivalent  (or  per  one 
heme  of  hemoglobin).  Such  confusion 
was  propagated  in  some  recent  publi¬ 
cations  [19],  including  our  own  studies 
[20]— [22],  where  the  s  values  given  by 
Zijlstra  et  al.  have  been  applied  [7]. 

It  gradually  becomes  clear  to  the 
authors  that  for  NIRS  of  hemoglobin 
quantification,  a  factor  of  four  needs 
to  be  multiplied  by  the  e  values  pub¬ 
lished  by  the  conventional  biochem¬ 
istry  methods.  This  is  a  conversion 
factor  to  account  for  four  hemes  per 
hemoglobin  molecule,  which  has  been 


During  the  last  ten  to 
15  years,  near  infrared 
spectroscopy  (NIRS) 
and  imaging  for  tissue 
vascular  oxygenation, 
breast  tumor  detection, 
and  functional 
brain  imaging  have 
been  intensively 
developed  for 
in  vivo  measurements. 


demonstrated  by  Mark  Cope  [10].  In 
this  way,  more  meaningful  results  for 
concentrations  of  [Hb]  and  [Hb02] 
can  be  arrived  from  the  NIR  absorp¬ 
tion  measurements. 


Furthermore,  accurate  determination 
of  [Hb]  and  [Hb02],  or  A[Hb]  and 
A[Hb02],  is  highly  sensitive  to  the 
accuracy  of  e  values,  and  intensive 
studies  have  been  conducted  over 
decades  on  how  to  determine  accurate 
e  values  [6] — [8],  [13] — [14].  The  e  val¬ 
ues  used  in  our  earlier  publications 
were  obtained  using  data  interpolation 
from  Zijlstra’s  e  table  [7].  However, 
such  data  interpolation  for  e  at  differ¬ 
ent  wavelengths  contained  potential 
errors  because  of  the  nonlinearity  of  s. 
Therefore,  a  calibration  procedure  was 
necessary  and  correction  factors  had  to 
be  applied  in  our  earlier  algorithm 
development  [21].  On  the  other  hand, 
Mark  Cope  has  generated  a  complete  e 
table  at  each  wavelength  between  650 
and  999  nm,  while  taking  four  hemes 
into  consideration.  Such  an  e  table  has 
been  recently  well  accepted  and  used 
[23]— [25]. 

Comparison  of  Published 
Extinction  Coefficients  for  NIRS 
Given  by  Two  Groups 

Since  Hoppe-Seyler  first  described  the 
visible  bands  of  oxyhemoglobin  in 
1862  [26],  numerous  reports  on  hemo- 


Fig.  1.  Absorption  spectra  of  oxy-  and  deoxyhemoglobin  from  Wray  (9), 
Cope  (10),  and  Zijlstra  (7),  (8). 
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Table  1 

Hemoglobin  extinction  coefficients  from  two  different  sources. 

Wave¬ 

length 

(nm) 

£Hb(mlVI 

'■cm-1) 

£Hbo2  (mM 

'■cm  ') 

Human  (7) 

Human  (8) 

Human  (10) 

Difference 

(Absolute) 

Human  (7) 

Human  (8) 

Human  (10) 

Difference 

(Absolute) 

750 

0.39 

0.37 

0.3865 

0.0165 

0.14 

0.15 

0.1374 

0.0126 

758 

0.354 

0.4205 

0.0665 

0.1596 

0.1494 

0.0102 

775 

0.29 

0.32 

0.3120 

0.008 

0.17 

0.18 

0.1759 

0.0041 

785 

0.276 

0.2494 

0.0266 

0.192 

0.192 

0 

800 

0.20 

0.21 

0.2099 

0.0001 

0.20 

0.21 

0.2163 

0.0063 

805 

0.20 

0.2037 

0.21 

0.2243 

840 

0.19 

0.1954 

0.25 

0.2768 

globin  absorption  spectra  have 
appeared  in  the  literature.  However, 
only  a  few  have  tabular  forms  of 
hemoglobin  extinction  coefficients.  In 
particular,  two  groups  have  reported 
tabular  forms  of  hemoglobin  extinction 
coefficients  in  a  broad  range  of  the 
NIR  spectrum.  Zijlstra  et  al.  reported  e 
values  for  human  blood  [7],  followed 
by  another  publication  for  updated 
human  e  values  with  a  comparison  to 
those  of  rat  blood  [8].  Moreover,  Cope 
has  generated  a  complete  £  table  at 
each  wavelength  between  650  and  999 
nm,  while  taking  four  hemes  into 
account  for  in  vivo  NIRS  measure¬ 
ments.  Cope  also  compared  his  own 
experiments  with  other  previous 
reports,  finding  that  the  isobestic  point 
of  hemoglobin  is  shifted  to  798  ±1.5 
nm  from  800  nm  reported  by  Horecker 
[1],  805  nm  by  Barlow  [27]  and  815 
nm  by  van  Assendelft  [28]. 

It  is  worthwhile  to  note  that  the 
determination  of  e  values  has  resulted 
in  a  range  of  accuracy  claims,  even  for 
values  reported  from  the  same  group. 
For  example,  according  to  Zijlstra  et 
al.,  a  value  of  eHb  at  560  nm  in  [29] 
was  12.72  (ndVR'cm-1),  while  it  was 
later  changed  to  13.09  (mM^'cm"1)  in 
[7];  a  value  of  £Hbo2  at  630  nm  was 
0.11  (mM^'cm-1)  as  given  in  [7], 
while  it  was  also  reported  as  0.14 
(mM_1cm_1)  in  [8].  From  Cope’s 
group,  their  published  reports  show  that 


£nb  and  £Hbo2  at  752  nm  are  1.574 
(mM^'cm-1)  and  0.608  (mM~1cm_1), 
respectively,  as  given  in  [9],  while  the 
corresponding  £  values  are  1.5938 
(mM-'cnr1)  and  0.5613  (mM_1cm_1) 
in  [10].  A  complete  comparison  of  e 
values  based  on  [7]— [  10]  by  Zijlstra  et 
al.,  Wray  et  al.,  and  Cope  is  plotted  in 
Figure  1,  showing  the  similarities  and 
differences  among  the  reported  extinc¬ 
tion  coefficients  of  oxy-  and  deoxyhe¬ 
moglobin.  Notice  that  the  values 
plotted  from  Zijlstra’ s  group  have  been 
multiplied  by  a  factor  of  four  to 
account  for  four  hemes  per  hemoglobin 
molecule  for  in  vivo  measurements.  It 
is  obvious  that  certain  variations 
among  different  sets  of  s  values  exist, 
while  the  basic  trends  are  the  same. 

Table  1,  on  the  other  hand,  shows 
the  tabular  form  of  hemoglobin  extinc¬ 
tion  coefficients  from  these  two 
groups.  Extinction  coefficients  in  the 
left  two  columns  of  £Hb  and  £Hbo2  are 
from  Zijlstra’ s  two  reports  [7],  [8]  and 
those  in  the  third  column  of  £Hb  and 
£Hbo2  are  from  Cope’s  report  [10]. 
Since  the  tabular  form  of  extinction 
coefficients  from  Zijlstra  does  not  have 
values  of  758  and  785  nm,  the  extinc¬ 
tion  coefficients  of  758  and  785  nm 
were  obtained  by  linear  interpolation 
between  750  and  775  nm  and  between 
775  and  800  nm,  respectively.  The 
hemoglobin  extinction  coefficients 
from  Cope  were  divided  by  four  to  be 


“per  equivalent”  so  that  they  could  be 
compared  to  those  from  Zijlstra  et  al. 
The  differences  listed  in  the  last  col¬ 
umn  are  the  differences  between 
Cope’s  £  values  and  the  values  from 
Zijlstra  et  al  [8]. 

While  the  differences  in  e  values 
seem  small  among  the  published  data, 
they  can  result  in  ~20%  variation  in 
quantification  of  hemoglobin  concen¬ 
tration,  i.e.,  [Hb],  [Hb02],  and 
[Hb]  totai  for  the  NIRS  determined  val¬ 
ues.  The  error  range  varies,  depending 
on  the  actual  wavelengths  used.  We 
are  currently  conducting  more  quanti¬ 
tative  error  analysis  to  demonstrate 
the  importance  of  accurate  e  for 
NIRS.  Meanwhile,  we  wish  in  this 
report  to  help  biomedical  engineers 
better  understand  the  extinction  coef¬ 
ficients  of  hemoglobin  and  to  further 
draw  attention  to  the  dependence  of 
NIRS  accuracy  on  the  extinction  coef¬ 
ficients  of  hemoglobin. 
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Abstract:  We  demonstrate  the  ability  to  investigate  breast  tumor  oxygen  dynamics  simultaneously 
by  steady-state  diffuse  re  flectance  s  pectroscopy  a  nd  F  OXY  oxygen  sens  or  i  n  response  t  o 
hyperbaric  oxygen  intervention. 

©2006  Optical  Society  of  America 

OCIS  codes:  (170.1470)  blood/tissue  constituent  monitoring;  (170.6510)  Spectroscopy,  tissue  diagnostics 

1.  Introduction 

Tumor  hy  poxia  has  pr  oven  t  o  have  a  prognostic  impact  in  cance  rs  a  nd  is  a  ssociated  w  ith  poor  r  esponse  to 
radiotherapy  [1,2].  Hyperbaric  o  xygenation  offers  ana  pproach  too  vercome  t  issue  ort  umor  hy  poxia.  Se  veral 
randomized  studies  investigating  hyperbaric  oxygenation  as  a  radiosensitizer  were  performed  30-40  years  ago  [3-5]. 
However,  hyperbaric  oxygenation  was  largely  abandoned.  Recently,  a  few  groups  have  developed  and  demonstrated 
a  new  st  rategy  to  combine  hyperbaric  oxygenation  and  radiotherapy,  based  on  the  finding  that  tumors  maintained 
elevated  oxygenation  even  after  exposure  to  hyperbaric  oxygen  [6-8].  Accordingly,  accurate  evaluation  of  tumor 
oxygenation  co  uld  be  era  cial  for  determination  of  radio-sensitivity.  In  our  st  udy,  we  ap  plied  two  m  onitoring 
techniques  simultaneously,  i  .e.,  steady-state  di  ffuse  reflectance  spectroscopy  (SSDRS)  to  monitor  turn  or  v  ascular 
oxygenation  and  a  m  ulti-channel  FOXY  oxygen  sensor  to  measure  regional  tumor  oxygen  tension  when  the  rats 
were  under  hyperbaric  oxygen  exposure.  The  objective  of  study  was  to  investigate  breast  tumor  oxygen  dynamics  in 
response  to  hyperbaric  ox  ygenation  and  t  o  exam  ine  wh  ether  a  p  ossible  im  provement  i  n  t  umor  o  xygenation  i  s 
maintained  after  hyperbaric  oxygenation. 

2.  Materials  and  Methods 

2.1  Animal  Preparation  and  experimental  setup 

Female  Fischer  344  rats  were  used  with  subcutaneously  growing  mammary  adenocarcinoma  13  762NF  (originally 
obtained  from  DCT,  NIH)  on  the  dorsum  of  the  thigh.  When  the  tumors  reached  ~1  cm  in  diam  eter,  the  rats  were 
anesthetized  with  ketamine  hy  drochloride  (  1.5ml;  100m  g/ml;  Aveco,  Fort  Dodge,  I  A)  a  nd  xy  lazine  by 
intraperitoneal  in  j  ection.  Turn  or  h  air  was  trimmed  fo  r  th  e  ease  of  optical  contact  fo  r  tran  smitting  NIR  ligh  t  an  d 
FOXY  probe  insertion.  The  rats  were  placed  in  the  hyperbaric  chamber  on  their  sides,  and  then  the  probes  of  SSDRS 
and  FOXY  were  fixed  securely  on  the  tumor.  Tumor  oxygenation  parameters  were  measured  simultaneously  by 
SSDRS  and  FOXY  during  respiratory  challenge  with  hyperbaric  oxygen  (HBO)  or  carbogen  (HBCB). 

During  the  experiments,  the  rats  were  exposed  to  hyperbaric  oxygen  in  an  acrylic  hyperbaric  chamber.  The  chamber 
was  flushed  with  air  for  15  min,  followed  by  normobaric  100%  oxygen  or  carbogen  for  15  min,  then  followed  by  a 
pressure  increase  to  2  atmospheres  absolute  (HBO)  for  30  min.  The  chamber  pressure  was  then  reduced  to  ambient, 
followed  by  a  flushing  with  air.  Gas  ex  change  (air  to  oxygen/carbogen)  was  accomplished  at  an  initial  flow  rate  of 
approximately  1  5  1/min.  Compression  and  decompression  ( from  1  t  o  2  atmosphere  oxygen/carbogen,  and  2  to  1 
atmosphere  oxygen/carbogen)  required  approximately  2  min. 

2.2  Steady-state  diffuse  reflectance  spectroscopy  (SSDRS)  for  measuring  changes  in  tumor  vascular 
oxygenation  (AHb02) 

The  instrument  used  to  acquire  reflectance  spectra  from  tumor  tissue  is  a  broadband  diffuse  reflectance  spectrometer. 
Briefly,  continuous  wave  (CW)  light  from  a  20  W  tungsten-halogen  light  source  (HL-2000HP,  ocean  optics,  FL)  is 
coupled  into  a  2.6-mm  core  diameter  fiber  optic  bundle,  the  distal  end  of  which  is  placed  in  physical  contact  with  the 
surface  of  the  tumor.  After  being  scatte  red  in  the  t  umor  tissue,  the  t  ransmitted  light  is  collected  by  a  1-mm  core 


1 


diameter  detection  fiber,  the  end  of  which  is  coupled  to  a  hand-held  spectrometer  (USB2000,  Ocean  optics,  FL).  The 
broadband  light  diffuse  spectrometer  provides  reflectance  spectra  from  400  to  900  nm. 

According  to  the  modified  Beer-Lambert  law,  as  given  in  Eqs.  (1)  and  (2),  changes  of  oxy-  and  deoxy-hemoglobin 
concentration,  A[Hb02]  and  A[Flb],  can  be  derived  from  the  measured  amplitudes  at  two  wavelengths  (750nm  and 
830nm),  by  using  extinction  coefficients  of  oxy-  and  deoxy-hemoglobin  published  by  Cope  [9]. 

logA™  =  (Aift750  A[ffi>]  +  SfibO-,  750  A[Tft02 1) '  DPF  ■  d  d) 

At 

logA830  =(£mmA[Hb]  +  £HboS30A[HbO2])DPF -d  (2) 

At 

where  Ah  is  the  baseline  amplitude,  A,  is  the  transient  amplitude  during  the  intervention,  and  d  is  the  direct  source- 
detector  separation.  DPF  (differential  path-length  factor)  is  a  tissue-dependent  parameter  and  defined  as  the  ratio 
between  the  optical  path  length  and  the  physical  separation  between  the  source  and  detector. 

2.3  FOXY  oxygen  sensor  for  measuring  oxygen  tension  of  tumors  (p02) 

Simultaneously,  regi  onal  p02  v  allies  in  tu  mors  were  monitored  by  a  m  ulti-channel,  fib  er  op  tic  ox  ygen  sen  sor 
(FOXY,  Ocean  Optics  Inc,  Dunedin,  FL).  Three  fluorescence-quenched,  optical  fiber  probes  (AL300,  tip  diameter 
410  pm)  were  i  nserted  into  different  regions  of  the  tumors.  Light  from  a  pulsed  blue  LED  (475  nm)  was  c  oupled 
into  one  branc  h  of  a  bifurcated  optical  fiber  probe  and  propagated  to  the  probe  tip.  The  distal  end  of  the  probe  is 
coated  with  a  th  in  layer  of  a  h  ydrophobic  so  1  g  el  m  aterial,  in  wh  ich  ox  ygen-sensing  ruthenium  co  mplex  i  s 
effectively  trapped.  Illumination  of  the  ruthenium  complex  causes  fluorescence  at  ~600  nm.  If  the  excited  ruthenium 
complex  encounters  oxygen  molecule,  the  excess  energy  is  transferred  to  the  oxygen  molecule,  thus  quenching  the 
fluorescence  signal.  The  degree  of  quenching  correlates  with  the  oxygen  concentration,  p02. 

3.  Results 
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Fig.  1  Time  course  of  A[Hb02]  and  pC>2  in  response  to  respiratory  challenge  for  a  representative  tumor,  simultaneous  measured  by 
SSDRS  (thicker  curve)  and  FOXY  (thinner  curves).  The  symbol  represents  the  time  needed  for  decompression. 


Figure  1  s hows  si multaneous  rea dings  o  f  c hanges  i n  [Hb02]  m easured  by  S  SDRS  and  three  p02  traces  in  three 
regions  of  a  breast  tumor  recorded  by  F  OXY  i  n  res  ponse  to  resp  iratory  eh  allenge.  This  t  umor  sho  wed  a  rapid 
response  to  oxyg  en  intervention  in  normal  pressure  with  a  significant  increase  in  A[Hb02].  The  three  readings  in 
p02  displayed  heterogeneous  responses  to  normobaric  oxygen  inhalation.  One  region  showed  distinct  improvement 
in  pO 2,  wh ile  th e  o ther  two  regions  had  little  o  r  no  improvement  in  pO 2.  Whe n  the  oxygen  was  pr  essured  to  2 
atmospheres,  A[Hb02]  had  a  further  increase  with  an  initial  increasing  rate  less  than  that  achieved  with  normobaric 
oxygen.  H  owever,  p  02  read  ings  in  th  e  three  resp  ective  reg  ions  sho  wed  sign  ificant  i  mprovement  un  der  th  e 
hyperbaric  oxygen  condition.  Returning  to  air  from  hyperbaric  oxygen  produced  a  gradual  decline  for  both  A[Hb02] 
and  pO  2-  A  strong  lin  ear  correlation  was  observed  between  maximal  Ap02  achi  eved  with  hy  perbaric 
Oxygen/Carbogen  intervention  and  that  achieved  at  the  same  location  with  respect  to  normobaric  oxygen/carbogen 
intervention  in  11  tumors  (Fig.  2). 
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Fig.  2  T  he  cor  relation  between  M  aximal  ApC>2 
(unit:  mmHg)  during  hyperbaric  Oxygen/C arbogen 
intervention  an  d  that  dur  ing  norm  obaric 
oxygen/carbogen  intervention  in  1 1  tumors. 

CB:  carbogen. 
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Fig.  3  Aver  age  p  O2  for  a  gr  oup  of  6  tu  mors  with 
hyperbaric  oxy  gen  inter  vention  achieved  with 
different  gas.  F1BO:  hyperbaric  oxygen. 


Tumor  pO  2  (a  verage  val  ue  of  a  group  of  si  x  t  umors  w  ith  HB  O)  i  n  a  ir  was  1 1 .7+2.6  mm  Hg,  and  it  increased 
significantly  d  uring  n  ormobaric  o  xygen  t  o  17.5+4  mm  Hg  (p  <0.0 1).  Hyperbaric  ox  ygenation  r  esulted  inam  ore 
significant  increase  in  pOi  than  normobaric  oxygenation  (p <0.005),  as  show  n  in  Figu  re  3.  Following  hyperbaric 
oxygen,  the  pCB  decrease  to  21.5+4.6  mmHg  in  10  min,  16.9+3.1  mmHg  in  15  min,  and  then  14.3+2.4  mmHg  in  20 
min.  Tumor  p02  values  showed  similar  responses  for  a  group  of  5  tumors  with  hyperbaric  carbogen  intervention. 


4.  Discussion 

This  study  sho  ws  that  th  e  two  techn  iques,  i.e.,  SSDRS  and  FOXY,  a  re  cons  istent  and  complementary  with  one 
another  fo r  t umor  o ximetry.  Our  data  demonstrates  that  hyperbaric  oxygenation  improves  the  tumor  oxygenation 
better  t  han  n  ormobaric  o  xygen.  It  al  so  su  ggests  t  hat  an  i  mprovement  in  t  umor  t  issue  oxy  genation  achieved  by 
hyperbaric  oxygenation  may  be  maintained  over  10-20  minutes  even  after  hyperbaric  oxygenation. 
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